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SUMMARY 


The  objective  of  this  project  was  to  establish  the  theoreti¬ 
cal  foundations  for  applying  higher  order  Rayleigh  and  interface 
modes  to  resonators  for  high  stability  sensors.  Existing  crystal 
resonators,  because  of  their  surface  sensitivity  and  thermal 
characteristics,  are  not  able  to  meet  the  stability  requirements 
of  Air  Force  tactical  systems. 

Theoretical  studies  of  the  charr*cteristics  of  Rayleigh  waves 
in  single  and  multi-layer  structures  were  carried  out  for  ST-cut 
Quartz  substrates.  Specified  materials  for  this  study  were  Mgo, 
Y  0  ,  AlN ,  SiO  ,  Al  0  ,  and  TiO  .  Velocity  dispersion  and  coupl- 
ing  to  bulk  modes  were  predicted  and  confirmed  by  comparison  with 
experimental  results  for  sputtered  films. 

Stoneley  ( interfacial )  waves  were  investigated  as  possible 
candidates  for  resonator  structures  because  of  their  inherently 
stable  and  low  cost  geometry  without  any  type  of  hermetic  enclos¬ 
ure  or  package.  A  successful  search  technique  was  developed  whi¬ 
ch  accurately  predicts  the  existence  of  Stoneley  waves  in  general 
anisotropic  (and  piezoelectric)  materials.  Using  this  technique, 
Stoneley  waves  were  for  the  first  time  predicted  in  single  crys¬ 
tal  quartz.  Several  useful  orientations  where  Stoneley  waves  are 
well  bound  and  piezoelectr ically  active  were  found  to  exist. 

Future  plans  are  to  design,  fabricate  and  test  interface 
wave  resonators.  The  end  goal  is  to  utilize  interface  wave  res¬ 
onators  and  sensors  in  Air  Force  tactical  systems. 


1.0  INTRODUCTION 


The  purpose  of  the  work  reported  was  to  perform  analyses  of 
the  propagation  characteristics  of  higher  order  acoustic  Rayleigh 
waves  in  multiple  layered,  anisotropic  materials.  The  objective 
was  to  establish  theoretical  foundations  for  predicting  the  per¬ 
formance  of  higher  order  acoustic  waves  and  to  apply  this  know¬ 
ledge  to  the  design  of  acoustic  resonators  as  strain  sensors  for 
cantilever  beam- type  accelerometers.  These  higher  order  Rayleigh 
waves  differ  from  the  normal  surface  waves  because  their  energy 
is  distributed  between  the  layers  and  the  substrate,  but  not  ent¬ 
irely  at  the  top  surface  of  the  layers.  This  has  the  advantage 
that  the  propagation  characteristics  are  less  sensitive  to  conta¬ 
mination  or  other  undesireable  effects  at  the  surface  of  the  cry¬ 
stal.  The  reflective  grooves  in  the  resonator  structure  can  be 
passivated  by  the  layers,  thus  reducing  the  possibility  of  false 

sensor  readings.  These  conditions  all  contribute  to  making  a  more 
stable  sensor  . 

1.1  Surface  Acoustic  Wave  Accelerometer 

A  Rayleigh  acoustic  wave  is  an  elastic  wave  which  propagates 
along  the  surface  of  an  elastic  solid.  Since  its  energy  is  con- 
confined  near  the  surface,  it  can  be  generated,  reflected,  and 
detected  on  the  surface,  Rayleigh  wave  devices  constitute  what 
is  called  a  planar  technology  and  hence  can  be  fabricated  by 
photolithographic  techniques.  Thus,  these  sensors  can  be  produc¬ 
ed  at  costs  comparable  to  those  of  silicon  integrated  circuits. 


The  Rayleigh  wave  resonator,  introduced  by  Dr.  Staples  in 
1974,  consists  of  a  pair  of  reflective  gratings  on  a  surface  as 
mirrors  in  a  Fabry-Perot  resonator  structure.  These  waves  can 
be  generated  and  detected  by  an  interdigital  transducer  on  a 
piezoelectric  substrate.  The  output  of  the  acoustically  resonant 
cavity,  shown  in  Figure.  1,  is  connected  to  a  feedback  amplifier 
to  maintian  a  high  Q,  standing  wave  in  the  cavity  formed  by  the 
two  grating  reflectors. 


Figure  1-  Typical  Rayleigh  wave  resonator  structure  showing  refl¬ 
ective  gratings  (acoustic  mirrors)  and  the  input/output 
interdigital  transducer. 

Since  Rayleigh  waves  are  confined  to  the  surface,  the  sensor 
can  be  configured  as  a  cantilever  beam,  as  shown  in  Figure  2.  In 
this  "push-pull"  open  loop  configuration,  thermal  effects  are  re¬ 
duced  by  using  two  resonators  (top  and  bottom)  and  the  output  is 


Figure 


2-  Cantilever  beam-type  rayleigh  Wave  Accelerometer  where 
the  output  is  the  difference  frequency  between  the 
opposing  top  and  bottom  resanant  cavities* 


1.2  Present  SAW  Technology  Limitations 

At  the  present  time,  short-term  frequency  stability  of  Ray¬ 
leigh  wave  resonators  is  comparable  to  that  of  bulk-wave  crystal 
oscillators.  Long-term  stability  of  fundamental  mode  Rayleigh 
resonators1 is  not  as  good  as  that  of  bulk  wave  crystal  oscilla¬ 
tors.  Rayleigh  wave  resonators  typically  have  short-term  (sec¬ 
onds  to  minutes)  frequency  stability  with  a  root  Allan  variance 

of  10"10  parts/part  of  the  resonant  freuqency,  and  long-term  ag- 

—  8 

•ing  rates  in  the  range  of  10  part/part/day . 

Short-term  stabilty  for  tactical  missile  applications  (fli- 
ght  time  typically  is  less  tan  20  minutes)  is  not  a  significant 
problem.  A  variation  of  m'!°  of  500  MHz  would  equal  0.05  Hz. 
With  an  assumed  quantization  of  0.1  fps/bit  and  input  accelera¬ 
tion  of  200  g,  this  would  equal  100  micro-g  acceleration  error. 
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Long-term  frequency  stability  is  a  more  serious  problem 
since  it  controls  the  bias  stability  of  Rayleigh  wave  sensors. 

The  dual  Rayleigh  wave  sensor  aging  rate  under  "ovenized"  condit- 

(2) 

ions  is  about  10-8  per  year.  For  a  Rayleigh  wave  accelerometer 
with  500  MHz  frequency  and  0.1  fps/bit,  this  gives  a  bias  stabil¬ 
ity  error  of  0.016  g/year.  This  stability  is  adequate  for  "moder¬ 
ately  accurate"  inertial  grade  accelerometer  applications.  With¬ 
out  oven  control  of  the  temperature,  the  oscillator  stability  is 
typically  10  6 /year  and  the  resulting  bias  stability  0.02  g/year 
which  exceeds  current  tactical  missile  specifications.  Tempera¬ 
ture  control  is  a  definite  limitation  of  current  Rayleigh  wave 
technology. 

An  important  limitation  of  Rayleigh  mode  resonators  is  that 

90%  of  its  energy  is  on  the  surface,  hence  it  is  very  sensitive 

to  surface  conditions.  Changes  of  1  part  per  billion  have  been 

(2) 

observed  m  the  resonator  frequency  in  a  period  of  one  hour.  This 
can  be  traced  to  molecular  surface  chemistry,  dust  particles,  or 
other  surface  effects  on  the  substrate.  After  the  SAW  pattern  is 
defined  it  would  be  highly  desireable  to  cover  or  passivate  the 
surface  so  as  to  protect  the  wave  from  surface  contaimination . 
This  would  negate  the  need  for  costly  packaging  and  sealing  of 
the  sensors.  This  last  limitation  is  addressed  by  the  current 
study  where  higher  order  acoustic  solutions  are  sought  which  do 
not  possess  the  high  degree  of  surface  sensitivity. 
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2.0  OBJECTIVES 


The  objective  of  the  phase  I  research  was  to  establish  the 
theoretical  foundations  for  using  higher  order  Rayleigh  and  int¬ 
erfacial  wave  resonators  as  strain  sensors  in  cantilever-type 
accelerometers.  Existing  crystals,  because  of  their  surface  sen¬ 
sitivity  and  thermal  characteristics,  are  not  able  to  meet  the 
technical  specifications  of  inertial  grade  accelerometers  for  Air 
Force  tactical  systems. 

The  overall  research  objective  is  to  obtain  resonator  struc¬ 
tures  with  improved  aging  characteristics,  so  that  these  sensors 
can  meet  the  performance  requirmenets  for  "moderate  accuracy" 
accelerometers  of  Air  Force  tactical  inertial  navigation  systems, 
including  missiles  and  aircraft. 

2.1  Work  statement 


The  research  c  ^ducted  under  Phase  I  was  defined  by  the 


following  work  statement: 


Task  1.  Perform  theoretical  studies  of  Rayleigh  waves  at  single 
layer  surfaces  on  quartz  substrates,  to  determine  the 
thickness  limitations  and  effects  of  the  layer.  Select¬ 
ed  materials  for  this  study  will  include  layered  acous¬ 
tic  wave  resonators  of  MgO,  Y20-r  AlN,  SiO  ,  Al203,  and 
Ti02. 

Task  2.  Correlate  the  theoretical  analysis  with  any  existing 
published  experimental  results  in  either  technical 
journals  or  DoD  reports. 


Task  3.  Expand  the  analysis  of  single  layers  to  three  anisotr¬ 
opic  layers.  Studies  will  be  made  in  characterizing 
various  combination  of  "loading  and  stiffening"  layers 
on  quartz. 

Task  4.  Design  and  develop  a  high  stability  Rayleigh  wave  sen — 
sor  using  the  theoretical  reults  of  the  above  studies. 


Task  5.  Investigate  theoretically  and  numerically  Stoneley  wav¬ 
es  in  joined  piezoelectric  half-spaces. 
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3.0  TECHNICAL  APPROACH 


Higner  order  Rayleigh  waves  differ  from  the  normal  surface 
waves  because  their  energy  is  distributed  between  the  layer  and 
the  substrate.  This  has  the  advantage  that  the  wave  characteris¬ 
tics  such  as  velocity  and  impedance  are  less  sensitive  to  changes 
in  the  surface  acoustic  and/or  electrical  conditions.  Also,  the 
reflective  grooves  and  electrode  materials  in  a  resonator  struc¬ 
ture  are  passivated  by  the  layer  and  the  effects  of  contamination 
are  reduced.  These  conditions  all  contribute  to  the  making  of  a 
more  stable  resonator. 

(3) 

Lord  Rayleigh  in  1887  predicted  the  propagation  of  seismic 

(A) 

waves  along  the  surface  of  the  earth.  Ewing,  Jardetzky  and  Press 

extended  Rayleigh's  work  to  layered  isotropic  media.  Since  most 

(5) 

crystals  are  anisotropic,  Lim  and  Musgrave  in  1970  studied  in¬ 
terface  wave  propagation  at  boundaries  between  two  similar  aniso- 

(6) 

tropic  media  and  they  extended  the  work  of  Stoneley  to  include 
anisotropic  media. 

(7) 

Farnell  and  Adler  in  1972  studied  Rayleigh  wave  propagation 
in  layered  anisotropic  media.  They  concentrated  their  studies  on 
special  cases  where  the  thickness  of  the  layers  was  small  compar¬ 
ed  to  a  wavelength.  The  results  showed  the  existence  of  single¬ 
mode  Rayleigh  waves  at  the  boundary  of  half  spaces  with  a  single 
layer  of  dissimilar  material. 

Recently  the  advantages  of  layered  wave  behavior  has  been 

(8) 

demonstrated  by  Yamonouchi  for  SiOx  on  LiTa(>3  and  LiNbC>3.  in 
addition  to  the  advantages  of  surface  passivation,  these  layered 
substrates  show  temperature  compensation  due  to  the  dissimilar 
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thermal  expansion  and  elastic  constants  with  temperature. 

In  order  to  study  propagation  of  waves  in  layered  acoustic 

media  two  approaches  were  taken.  The  first  followed  the  analysis 
(9) 

of  Farnell  to  describe  and  predict  wave  behavior  for  Rayleigh 
wave  propagation  in  layered  media.  This  analysis  was  used  to 
evaluate  several  layer  materials  deposited  on  ST-cut  quartz.  The 
theoretical  predictions  were  then  compared  with  previously  publ¬ 
ished  experimental  results. 

The  second  approach  was  undertaken  jointly  with  Professor 
D.  Barnett  of  Stanford  University.  Professor  Barnett  specifi¬ 
cally  analyzed  search  methods  for  interfacial  wave  propagation 

in  general  anisotropic  media.  These  search  methods  were  derived 

(10) 

from  the  theory  published  by  Barnett  and  Lothe.  A  reprint  of 
this  paper  is  included  with  this  report  as  Appendix  II.  Theory 
of  interfacial  wave  propagation  is  described  in  section  3,  and 
numerical  results  for  prescribed  materials  are  presented  in  sect¬ 
ion  4.  The  importance  of  this  work  is  that  a  method  was  develop¬ 
ed  which  predicts  the  existence  of  interface  waves.  Unlike  sur¬ 
face  waves,  interface  waves  do  not  exist  for  all  material  combin¬ 
ations.  Until  now  an  exhausting  search  was  always  required  to 
determine  the  existence  of  such  waves.  Using  new  search  techni¬ 
ques,  Stoneley  wave  behavior  has  been  predicted  without  complet¬ 
ely  solving  wave  equations.  For  the  first  time  Stoneley  waves 
have  been  found  in  quartz.  By  properly  selecting  the  cuts  of 
quartz  used,  it  was  found  that  Stoneley  waves  with  good  piezo¬ 
electric  coupling  and  fast  decay  with  distance  from  the  interface 
do  exist. 
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3.1  Surface  Acoustic  Waves 


In  a  piezoelectric  material  the  three  displacement  components 
Uj  along  the  Cartesian  axes  x..  ,  and  t’  electrostatic  potential 
are  related  by  the  equations  of  motion  and  electrostatics. 


d2Ui  _  d2uk 

P  H  2  Ci)U  ~ 

at  dXi  dXi 


djp 

dXi  dxk 


=  0 


d2Uk  d2<p 

dXi  dxt  dXi  dxk 


(1) 


in  which  p'  is  the  density  of  the  material  and  c,  ,  ,  e  and 

ijkl  kij 

are  respectively,  the  tensors  for  the  elastic  stiffness  mea¬ 
sured  at  constant  electric  field,  for  the  piezoelectric  coupling 
and  for  dielectric  permittivity  measured  at  constant  strain.  The 
subscripts  i,j,k,  and  1  range  over  the  values  1,2,  and  3  anpl  the 
summation  convention  on  repeated  subscripts  is  employed.  The 
contracted  matrix  notation  and  the  full  tensor  notation  will  be 
used  interchangeably  in  this  report.  Particular  elastic,  piezo¬ 
electric,  and  dielectric  constants  for  materials  studied  in  this 
report  are  given  in  Appendix  I. 

All  solutions  considered  are  straight  crested  in  that  the 
planes  of  constant  phase  are  perpendicular  to  the  direction  of 
propagation  and  there  is  no  variation  of  the  amplitude  in  the 
direction  perpendicular  to  the  sagittal  plane.  The  coordinate 
system  used  throughout  this  report  is  as  shown  in  Figure  3  with 
the  direction  taken  as  parallel  to  the  propagation  vector  and 

the  X3  direction  taken  as  perpendicular  to  the  surface  (also  in¬ 
terface  in  later  sections) . 

/ 


9 


*3 


Figure  3-  Diagram  showing  coordinate  system  orientation  for 
theoretical  analysis  of  wave  propagation. 

For  a  single  material  surface,  Rayleigh  wave  solutions  will 
be  of  the  form 


Uj  =  cij  exp  (ikbx*)  exp  [ik(x:  —  vi)] 
<p  —  di  exp  (ikbx-i)  exp  [ik(x.i  —  vt)]. 


Substituting  (2)  into  the  equations  of  motion  (1)  results  in 
secular  matrix  of  the  form 


A  13 

r,« 


rM 

Tu 


(Tn  pv  )  r  12  r*i3 

Ti2  (1^22  “  PV  )  r*23 

"  r23  (r33  -  p*-’2) 

r24  r34.  r44 


•  ' 

<*1 

a2 

•«4< 

=  o. 


(3) 


where  the  terms  of  the  matrix  are  related  to  the  material  const 
ants  by  the  following  quadratic  coefficients: 


ra  =  4"  2ciu36  c»ui 

Py*  =  £2323^  ”1”  2Cjaa3i»  4“ 

^33  ==  ^3333^  4*  -^1333^  4”  ^13J3 

Pl2  888  Cj323 &  4*  1123  4-  Ct3u)k  4“  Cn)2 

rl3  =  cu33k2  4-  (cn 33  4-  cl3 i3)6  4*  Cjju 

P23  —  C3333^  4*  (^1233  4“  ^23 13)^  4"  Cl213) 

Tli  —  *331^  4-  (e  131  4-  £311)^  4-  eux 

Ta*  =  C332i»2  4~  (fim  4-  ^312)^  4“  £ua 

r34  =  C33362  4~  (^133  4“  ^313 )k  4*  ^113) 

ri4  =  ~(«33^  4"  2cu^  4-  <ll)- 


The  above  homogeneous  equations  have  a  nontrivial  solution  only 
if  the  determinant  of  the  square  T  matrix  is  zero.  A  computer 
program  was  written  to  expand  the  above  matrix  assuming  a  known 
velocity,  v,  as  a  polynomial  (8th  order  for  piezoelectric  mater¬ 
ials)  in  the  decay  constant  b.  The  program  evaluated  the  complex 
roots  b±  which  occured  as  conjugate  pairs.  In  the  substrate  the 
solutions  of  interest  decay  with  depth  away  from  the  surface, 
hence  only  roots  in  the  lower  half  of  the  complex  plane  are  re¬ 
tained  . 


Im  aw  v*  y\ ».»  cjLi  rj^u  LV»o  ■£  v  o  O  c*nv  t1  OHS 

ii  Ui-UCi  lO  uww  batioty  uu'-  j-  a  ^  Du*  i.uv^  ^  jr  - - - ’*•“ 


a  wave  solution  composed  of  sums  of  terms  given  by  eq. (2)  where 
each  term  satisfies  the  wave  equation  (1)  is  used.  Each  term  is 

(n)  ,  . 

considered  an  eigenfunction , a  ,  with  an  associated  eigenvalue 
or  decay  constant  bM  .  The  total  solution  is  then  given  by 
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(5) 


ui  ~  12  '0% a\n)  exp  ( ikb(n)Xz )  exp  [ik(xx  —  vl )] 

V  ~  12  Ona\n)  exp  (ikb{n)x-4)  exp  [ik(xx  —  y£)] 

where  summation  over  all  the  root  numbers  (4  for  piezoelectrics)  . 

The  complete  solution  is  obtained  by  satisfying  the  mechani¬ 
cal  and  electrical  boundary  conditions.  Mechanically  the  normal 

stress  compoents  T  /  T  ,  and  T  must  be  zero  at  the  surface. 

13  23  23 

Electrically  there  must  be  continutity  of  electric  potential  and 
the  normal  componentf  displaceme.it  current  D  ^ 

Substitution  of  eqs.(5)  results  in  a  different  set  of  four 
homogeneous  equations  in  the  four  weighting  functions  Cn  ,  and 
for  nontrivial  solutions  the  determinant  of  the  coefficients  of  C 
must  vanish. 


j  (c«»i  4~  Cni*bH)a%n  4~  (eua  4*  <tyubn)cun 

(Cjuii  4~  +  (^ui  +  Urubn)txiu 

4"  4"  (tf i32  4~  tyub*)ccin 

“ '(£'#»»  +  CiiJ)n)atih  4*  ( 4*  tub”  lt.))aun 

L 
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An  objective  of  this  study  was  to  determine  the  behavior  of 
Rayleigh  waves  on  ST-cut  quartz.  Equations  1-6  were  solved  num¬ 
erically  using  an  iterative  search  on  velocity  until  equation  (6) 
was  satisfied.  A  plot  of  the  displacement  components  and  potent¬ 
ial  are  shown  in  Figure  4a.  When  the  surface  is  shorted  the  last 
row  of  eq. (6)  is  replaced  by  the  the  eigenvectors  which  forces 
the  potential  to  zero  at  the  surface.  A  plot  of  the  displacements 
and  potential  for  this  case  is  shown  in  Figure  4b. 

3.2  Surface  Acoustic  Waves  in  Layered  Media 

When  the  surface  of  a  crystal  is  covered  with  a  layer  of 
dissimilar  material  as  shown  in  Figure  5,  the  behavior  of  the 
Rayleigh  wave  is  modified  to  take  into  account  (1)  the  elastic 
properties  of  the  layer  and  (2)  the  electric  boundary  conditions 
of  the  layer  surface  and  the  interface  region. 

Solutions  to  the  layered  problem  can  be  obtained  by  solving 

equations  1-3  in  the  layer  and  substrate  separately  using  the 

appropriate  material  constants.  The  individual  solutions  in  the 
layer  and  infinite  substrate  are  still  assumed  to  be  the  same 

form  as  equation  5  where  the  appropriate  layer  or  substrate  roots 
apply. 

The  complete  solution  is  obtained  by  satisfying  the  approp¬ 
riate  boundary  conditions.  In  the  case  of  layered  media  the 
boundary  condition  determinant  is  more  complex.  The  surface  of 
the  layer  is  stress  free  and  electrically  open  or  shorted.  How¬ 
ever,  the  interface  region  must  also  be  satisfied  by  matching 
the  normal  stresses,  the  displacements,  and  electrical  potential 
and/or  displacement  currents. 
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NORMALIZED  DISPLACEMENT /POTENTIAL  NORMALIZED  DISPLRCEMENT/POTENTIAL 


X3 


Figure  5  -  Diagram  showing  coordinate  system  orientation  for 
theoretical  analysis  of  wave  propagation  in 
layered  media. 

In  this  report  three  distinct  types  of  materials  were  anal¬ 
yzed.  The  first  situation  involved  a  piezoelectric  substrate 

(YX-Quartz)  and  a  non-piezoelectric  isotropic  layer  such  as  fused 

(9 

quartz,  hereafter  designated  PN53  after  the  notation  of  Farnell. 
In  this  case  the  transverse  displacement  component  is  uncoupled 
and  there  are  7  boundary  conditions  which  must  be  satisfied.  The 
boundary  condition  equations  for  a  unshorted  interface  are 


T3(X3"0) 

1 

-t|<X3»0) 

T5<x3“0> 

1 

-T§(X3-0> 

T3<X3«h) 

I 

0 

T5<X3«h> 

i 

i 

0 

u^(X3-0) 

i 

1 

-uf <X3-0) 

u3CX3«0) 

l 

— \ 

1 

l 

-uf(X3-0) 

D^cxa-0) 

0 

(7) 
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and  for  a  shorted  interface  are 


t^cx3-0) 

-t|(X3-0) 

T5<X3-0> 

-t|cx3-0) 

T3<X3-h) 

0 

T5<X3-h> 

0 

U^Xg-O) 

-uf (X3«0) 

u3(X3-0) 

-u|<X3-0) 

♦LCX3-0)  I  0 


=  0 


(8) 


The  displacement  vectors  of  equations  7  and  8  are  modified  to 
account  for  the  dielectric  constant  of  the  layer  and  the  free 
space  above  the  layer.  When  the  substrate  crystal  solution  con¬ 
tains  transverse  displacements  these  must  be  accounted  for  in  the 
boundary  conditions.  This  applies  whether  the  layer  material  is 
isotropic  (no  coupled  transverse  displacements)  or  anisotropic. 
This  case  is  represented  by  ST-Quartz  substrates  and  layer  mater¬ 
ials  such  as  fused  quartz.  Rutile,  Sapphire,  and  yittrium  Oxide. 
In  these  cases  the  unshorted  boundary  conditions  are 


T3<X3-0) 

1 

-t|<x3-0) 

T5<X3-0) 

1 

-t|<x3-0) 

T^(X3-0) 

1 

-T*<X3“0) 

T3<X3-h) 

1 

I 

0 

T5<X3-h) 

1 

0 

T^(X3«h) 

• 

1 

0 

ui<X3-0) 

! 

-u®CX3-0> 

u3(X3=»0) 

1 

1 

-u|(X3=0) 

u3<X3-0) 

1 

-u|(X3-0) 

D3(X3»0) 

"  1  ’ 
I 

0 

=  0 


(9) 
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and  for  shorted  electrical  conditions  are 


T3<X3-0> 

Tg  <X3“0) 

1 

1 

-Tf(X3-0) 

-TfcXg-0) 

T4<X3-0) 

1 

-T*(X3-0) 

T^CX3«h) 

1 

1 

0 

T5<X3-h) 

1 

0 

T4<X3-h> 

1 

0 

u'fCXa-O) 

1 

-uf CX3-0) 

u3<X3-0) 

1 

| 

-uf (X3-0) 

u3(X3-0> 

1 

-u3(X3-0) 

4>LCX3«0) 

'  1 
! 

0 

T1 


=  0 


(10) 


The  boundary  conditions  are  further  complicated  when  the 


layer  material  is  itself  piezoelectric  and  with  coupled  tran¬ 
sverse  displacements.  This  case  was  analyzed  for  the  substrate 
material  ST-Quartz  and  layer  Aluminum  Nitride.  The  boundary 


condition  determinant  in  this  case  contains  twelve  equations. 


The  form  for  unshorted  surface  and  interface  is 


t^<x3»0) 

1 

-T3(X3*=0) 

t^(x3*0) 

1 

-T5<X3~0> 

T4<X3“0) 

1 

-tJ(X3»0> 

T3(X3=h) 

! 

| 

0 

T5<X3-h> 

1 

0 

T^CXg-h) 

1 

0 

u*J’(X3=0) 

1 

-uf (X3“0) 

u3(X3»0) 

1 

i 

-u|(X3«0) 

u3(X3-0) 

1 

-uf (X3-0) 

D3(X3=h) 

"  1  • 
1 

0 

/•(X3=0) 

1 

-«»S(X3=>0) 

D3(X3-0) 

! 

-Df (X3-0) 

=  0 


(11) 
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When  the  top  surface  and  interface  regions  are  shorted 


electrically  the  new  for  of  the  boundary  matrix  is 


T3<X3-0> 

1 

-t|(X3-0) 

Ts<X3-0) 

1 

-t|cx3-0) 

t1^(X3-0) 

1 

-T*<X3-0> 

T3<X3-h> 

1 

| 

0 

Ts<X3-h> 

1 

0 

T4<X3-h> 

1 

0 

Ul(X3«*0> 

1 

-uf CX3-0) 

u3(X3-0) 

1 

I 

-uf(X3-0> 

u3<X3**0) 

1 

-u3CX3-0> 

♦3<X3-h> 

“  1 

1 

0 

«L(X3-0) 

1 

0 

0 

1 

♦f(x3-0) 

In  the  general  case  there  will  be  a  solution  which  satisfies 
the  wave  equations  in  the  layer  and  substrate  and  simultaneously 
satisfies  all  12  of  the  above  boundary  conditions.  Electrical 
coupling  is  estimated  by  comparing  the  shorted  and  unshorted  wave 
velocities.  For  a  piezoelectric  layer  on  a  piezoelectric  subst¬ 
rate  there  will  be  four  possible  combinations  of  boundary  condit¬ 
ions  corresponding  to  four  possible  transducer  electrode  .‘onfig- 
urations  shown  in  Figure  6. 


_ oja_n-Ei _ 

fi.rf  r»  «  _ 

n  r«  n  n 

Figure  6-  Four  possible  transducer  electrode  configurations 
when  both  layer  and  substrate  are  piezoelectric. 
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3.3  Interfacial  Waves 


Interfacial  waves,  commonly  called  Stoneley  Waves,  occur  at 
the  boundary  of  two  infinite  dissimilar  materials  as  shown  in 
Figure  7.  Unlike  surface  acoustic  waves,  Stoneley  waves  only 
exist  under  special  circumstances.  The  exact  nature  of  these 
special  circumstances  is  not  well  understood  and  an  objective  of 
this  study  was  to  develop  methods  of  searching  for  Stoneley  waves 
in  general  anisotropic  media. 


*3 


Figure  7-  Diagram  showing  coordinate  system  orientation 
for  analysis  of  wave  propagation  along  the 
interface  between  two  anisotropic  materials. 


The  analysis  of  interfacial  waves  is  similar  to  analyzing 
layered  surface  acoustic  waves.  The  equations  of  motion  are 
solved  independently  in  each  material  and  then  summation  solut¬ 
ions  are  used  to  satisfy  the  boundary  conditions.  For  interface 
waves  the  unshorted  boundary  condition  matrix  is 
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T3(X3-0) 

1 

-T3<X3-0> 

T5<X3«0> 

1 
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OtCX3-0) 

“T 

t 

-*b(X3-0> 

D3CX3-0) 

l 

-D3(X3-0) 

=0 


(13) 


The  above  notation  refers  to  the  top  (t)  and  bottom  (b)  mat 
erials  and  their  respective  material  constants.  When  the  inter¬ 
face  is  shorted  by  a  conducting  plane  the  boundary  condition  mat 
rix  takes  the  form 


T3(X3“0) 

j  -T3(X3-0) 

T^(X3-0) 

1  -T^CX3-0) 

tJ(X3-0) 
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- - f  - 
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=0 


(14) 
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3.4  Multi-Layered  Resonators 


Multi-layered  resonators  are  photo  etched  on  crystals  cont¬ 
aining  multiple  layers  of  dissimilar  materials  on  a  piezoelectric 
substrate.  The  layers  may  be  piezoelectrically  active  also  and 
shorting  planes  may  exist  at  any  interface  as  depicted  in  figure 
8  for  a  two  layer  system. 


*3 


Figure  8-  Coordinate  system  for  a  two-layered  resonator 
geometry. 

Analysis  reveals  the  solution  contains  both  Rayleigh  and 
Stoneley  wave  like  solutions.  Because  of  this,  the  behavior  of 
Rayleigh  and  Stoneley  waves  in  the  individual  material  sub-sets 
was  studied  before  solutions  in  multi-layered  media  were  under¬ 
taken  . 

The  analysis  follows  a  familiar  trend.  The  equations  of 
motion  are  solved  independently  in  each  material  and  summation 
solutions  are  used  to  satisfy  the  boundary  conditions.  For  waves 


in  unshorted  multi-layered  resonators  the  boundary  condition  mat¬ 
rix  is  as  follows: 
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(15) 


The  above  notation  refers  to  the  top  layer (m) ,  intermediate 

layer (n)  and  substrate  (s) .  When  the  two  interfaces  and  top  layer 
are  shorted  the  boundary  condition  matrix  takes  the  form: 


22 


T3CX3-0) 

1 

-t§<x3«0> 

1  0 

T5<X3-0> 

1 

-T§(Xg-0) 

1  0 

t5<X3-0) 

1 

-t5cx3-0) 

1  0 

u*(X3-0) 

1 

-u" (X3-0) 

1  0 

U2<X3«0) 

1 

1 

-U2<X3-0) 

5  0 
| 

u|(X3-0) 

1 

| 

-U3CX3-0) 

;  0 

0 

1 

T3  ^X3"h 1 ) 

,-t5(x3-  hi) 

0 

1 

Tg  ^X3*h 1 ) 

|-Tg(X3-  hi) 

0 

1 

T^Xg-hl  ) 

I-tJ(X3-  hi) 

0 

0 

1 

1 

uf CX3-hl ) 
U2<X3«hn 

l-ulJ'(X3-  hi) 
'-UgCXg-  hi) 

0 

1 

1 

u^Xg-hl ) 

J-U3CX3-  hi) 

0 

1 

0 

1  T^(X3-  h2 ) 

0 

1 

0 

1  T^CX3-  h2 ) 

0 

1 

- 1- 

0 

1  tJ(X3“  h2 ) 

H - 

OS(X3«0) 

1 

0 

I  0 

0 

1 

*n<X3-0) 

1  0 

0 

1 

♦n(X3»hl) 

1  0 

0 

1 

1 

0 

*mCX3“  hi) 

0 

\ 

\ 

0 

,  <J>mCX3-  h2 ) 

Clearly  intuitive  methods  of  solution  determination  would  be 
difficult  in  this  case  due  to  20  equations  in  20  unknowns.  How¬ 
ever  recognizing  certain  groupings  as  natural  Rayleigh  and/or 
Stoneley  wave  solution  sub-sets  may  be  used  as  a  search  method 
in  these  multi-layered  structures. 
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4.0  TECHNICAL  RESULTS 


In  this  section  the  technical  results  of  research  conducted 
under  Phase  I  (Section  2.1)  is  presented.  In  Section  4.1  the 

results  of  layered  media  surface  wave  velocity  studies  are  pre¬ 
sented  for  six  different  layer  materials.  In  Section  4.2  the 
results  of  studies  on  interfacial  waves  are  presented.  The  dis¬ 
covery  of  interface  waves  (Stoneley)  in  quartz  for  the  first  time 
is  an  important  result  of  the  phase  I  study  and  may  lead  to  a  new 
class  of  acoustic  resonators.  Using  these  results,  a  SAW  sensor 
design  is  discussed  in  section  4.3  and  in  Section  4.4  theoretical 
results  are  compared  to  published  experimental  results  for  layer¬ 
ed  SAW  resonators. 

4.1  Layered  Substrate  Analyses 

In  this  section  six  different  materials  are  considered  as 

candidates  for  passivation  of  SAW  resonators.  The  materials 
considered  are  magnesium  oxide  (MgO) ,  yittrium  oxide  (Y2  0^  ), 

silicon  dioxide  (Si02) ,  aluminum  nitride  ( AIN) ,  rutile  or  titan¬ 
ium  dioxide  (Ti02  ) ,  and  sapphire  (Al2°3  ^  *  The  t^eoret^cal  aPPr“ 
oach  for  analyzing  layered  SAW  velocity  was  discussed  in  Section 
3.2.  Results  are  presented  as  dispersion  diagrams,  surface  wave 
coupling  (Av/Vo)  and  displacement/potential  vectors  vs  distance 
from  the  top  of  the  layered  surface.  For  completeness  tables  of 
numerical  solutions  are  included  with  displacement  diagrams. 
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4.1.1  Analysis  of  MgO/ST-Quartz 
Layers  of  MgO  deposited  on  ST  cut  Quartz  substrates  were 
analyzed.  Elastic  and  dielectric  constant  data  for  single  crys¬ 
tal  MgO  (Appendix  I)  were  used  with  the  crystal  c-axis  perpendic¬ 
ular  to  the  surface  and  the  X-axis  parallel  to  the  Quartz  X-axis. 

Surface  wave  velocity  vs  normalized  layer  thickness  is  shown 
in  Figure  9a.  Because  MgO  is  a  faster  (higher  velocity)  material 
the  velocity  increases  with  layer  thickness  until  reaching  the 
slow-shear  wave  velocity  (dotted  line)  of  ST-Quartz.  Surface  wave 
coupling  is  proportional  to  Av/Vo  and  this  parameter  is  plotted 
in  Figure  9b  for  two  transducer  configurations.  Coupling  drops 
quickly  for  layers  greater  than  1%  of  a  wavelength. 

Selected  numerical  solutions  for  MgO  layers  are  shown  in 
Figures  10,11,12,  and  13.  Figures  10-12  show  the  effect  of  short¬ 
ing  the  electric  field  when  the  layer  is  10%  of  the  wavelength. 
The  displacements/potential  show  penetration  increasing  with 
layer  thickness.  In  Figure  13  displacements  and  potential  are 
shown  for  a  layer  only  0.1%  of  the  wavelength.  The  increased 
penetration  with  increasing  layer  thickness  is  due  to  the  near¬ 
ness  of  the  bulk  shear  wave  in  the  Quartz  substrate. 

Numerically  MgO/ST-Quartz  is  governed  by  equations  9  and  10 
of  section  3.2  which  contain  10  sub-solutions.  These  are  present¬ 
ed  in  the  Tables  of  Figures  10-13.  Displacement/potential  vectors 
associated  with  each  sub— solution  root  (decay  constant)  are  shown 
together  with  the  10  complex  weighting  factors  for  the  complete 
summation  solution  per  equation  5  of  Section  3.1. 
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Figure  11-  Layered  velocity  solution  for  MgO/ST-Quartz 

when  the  top  layer  surface  is  shorted.  Layer 
thickness  is  10%  of  a  wavelength. 
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Figure  12-  Layered  velocity  solution  for  MgO/ST-Quartz 

when  the  top  layer  and  interface  are  shorted. 
Layer  thickness  is  10 %  of  a  wavelength. 
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Figure  13-  Layered  velocity  solution  for  MgO/ST-Quartz 

when  the  top  layer  and  interface  are  unshorted. 
Layer  thickness  is  0.1%  of  a  wavelength. 


4.1.2  Analysis  of  Y203  /ST-Quartz 

Layers  of  yittrium  oxide  on  ST-Quartz  substrates  were  analy¬ 
zed.  Elastic  and  dielectric  constant  data  for  single  crystal  ^O;} 
(Appendix  I)  were  used  with  the  crystal  c-axis  perpendicular  to 
the  surface  and  the  X-axis  parallel  to  the  Quartz  X-axis. 

Surface  wave  velocity  vs  normalized  layer  thickness  is  shown 
in  Figure  14a.  Because  Y203  has  a  velocity  only  slightly  high¬ 
er  than  that  of  the  quartz  substrate ,  the  velocity  of  the  layered 
wave  becomes  constant  for  large  normalized  thickness.  This  impli¬ 
es  the  wave  is  entirely  in  the  film  and  not  affected  by  the  subs¬ 
trate.  This  is  confirmed  by  the  coupling  shown  in  Figure  14b 
where  the  coupling  becomes  zero  as  the  thickness  approaches  a 
wavelength. 

Selected  numerical  solutions  for  Y203  layers  10%  of  a  wave¬ 
length  are  shown  in  Figures  15-17  under  different  conditions  of 
shorted  and  unshorted  electric  fields  at  the  top  layer  and  inter¬ 
face.  Numerically  Y2C>3  /ST-Quartz  is  governed  by  equations  9  and 
10  of  section  3.2  which  contain  10  sub-solutions.  These  are  pre¬ 
sented  in  the  Tables  of  Figures  15-17.  Displacement/potential 
vectors  associated  with  each  sub-solution  root  (decay  constant) 
are  shown  together  with  the  10  complex  weighting  factors  for  the 
complete  summation  solution  per  equation  5  of  section  3.1. 
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Figure  15-  Layered  velocity  solution  for  Y203/ST-Quartz 

when  the  top  layer  and  interface  are  unshorted. 
Layer  thickness  is  10%  of  a  wavelength. 


Figure  16-  Layered  velocity  solution  for  3^T-^uartz 
when  the  top  layer  surface  is  snorted.  Layer 
thickness  is  10%  of  a  wavelength. 
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4.1.3  Analysis  of  SiOx/ST-Quartz 

Layers  of  fused  quartz  on  ST-cut  Quartz  substrates  were  an¬ 
alyzed.  Elastic  and  dielectric  constant  data  for  isotropic  fused 
quartz  (Appendix  I)  were  used.  Surface  wave  velocity  vs  normal¬ 
ized  layer  thickness  is  shown  in  Figure  18a.  Because  fused  quartz 
is  a  faster  material  the  velocity  increases  with  layer  thickness 
until  reaching  the  slow  shear  velocity  of  the  substrate.  Surface 
wave  coupling  is  plotted  in  Figure  18b  for  two  transducer  config¬ 
urations.  The  coupling  becomes  zero  for  thicknesses  greater  than 
half  a  wavelength. 

Selected  numerical  solutions  for  fused  quartz  layers  are 
shown  in  Figures, 19-21 .  The  displacements/potential  show  penet¬ 
ration  increasing  with  layer  thickness.  The  increased  penetrat¬ 
ion  with  increasing  layer  thickness  is  due  to  the  nearness  of  the 
bulk  shear  wave  in  the  Quartz  substrate. 

Numerically  SiOx/ST-Quartz  is  governed  by  equations  9  and  10 
of  section  3.2  which  contain  10  sub-solutions.  Displacement  and 
potential  vectors  associated  with  each  sub-solution  are  shown  to¬ 
gether  with  the  10  complex  weighting  factors  for  the  complete 
summation  solution  per  equation  5  of  Section  3.1. 
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4.1.4  Analysis  of  AlN/ST-Quar tz 

Layers  of  aluminum  nitride  ( AIN)  deposited  on  ST-cut  Quartz 
substrates  were  analyzed.  Surface  wave  velocity  vs  normalized  la¬ 
yer  thickness  is  shown  in  Figure  22a  and  23a  for  combinations  of 
shorted  layer  and/or  interface.  Because  AIN  is  a  faster  material, 
the  velocity  increases  with  layer  thickness  until  reaching  the 
slow-shear  velocity  (dotted  line)  of  the  ST-Quartz  substrate.  Be¬ 
cause  AlN  is  piezoelectric,  four  possible  transducer  configurat¬ 
ions  are  possible.  Coupling  efficiency  to  surface  waves  with  only 
interdigital  electrodes  at  either  the  interface  or  top  layer  sur¬ 
face  is  shown  in  Figure  22b.  When  the  interdigital  electrodes  are 
at  the  interface,  the  coupling  is  reduced  because  the  wave  displ¬ 
acements  are  located  at  the  top  surface.  If  the  electrodes  are  at 
the  surface  and  the  layer  is  thick,  coupling  efficiency  improves 
because  of  the  piezoelectric  AIN  layer. 

For  the  floating  ground  plane  (either  at  the  surface  or  int¬ 
erface)  the  coupling  is  shown  in  Figure  23b.  Normally,  the  coup¬ 
ling  is  low  for  thin  layers.  However  when  the  layer  thickness  is 
increased,  the  coupling  peaks  at  2%  of  a  wavelength. 

Numerical  solutions  for  AlN  layers  are  shown  in  Figures  24- 
27  for  different  shorting  conditions  and  thickness  equal  to  10% 
of  a  wavelength.  Numerically  AlN/ST-Quar tz  is  governed  by  equat¬ 
ions  11  and  12  of  section  3.2  which  contain  12  sub-solutions. 
These  are  presented  in  the  Tables  of  Figures  24-27.  Displacement 
and  potential  vectors  associated  with  each  sub-solution  root  are 
shown  together  with  the  12  complex  weighting  factors  for  the  com¬ 
plete  summation  solution  per  equation  5  of  section  3.1. 
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DELTA  V/V  SRW  VELOCITY  (M/SEC) 


DELTA  V/V  SAW  VELOCITY  (M/SEC) 


SUBSTRATE: ST-CUT  QURRTZ 
LAYER: ZX  RL  NITRIDE  ■ 


NORMALIZED  LAYER  THICKNESS  (H/Lambda> 


SUBSTRATE: ST-CUT  QUARTZ 
LAYER: ZX  AL  NITRIDE 
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0. 0000 
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NORMALIZED  LAYER  THICKNESS  (H/Lambda) 

(b) 

Figure  23-  Velocity  dispersion  (a)  and  wave  coupling  (b) 
as  a  function  of  layer  thickness  for  AIN  on 
ST-Quartz  with  interdigital  electrodes  and  an 
opposing  ground  plane. 
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4.1.5  Analysis  of  Ti02/ST-Quartz 

Layers  of  rutile  (Ti02)  on  ST-cut  Quartz  substrates  were  an¬ 
alyzed.  Elastic  and  dielectric  constant  data  for  rutile  were  us¬ 
ed  (Appendix  I) .  Surface  wave  velocity  vs  normalized  layer  thick¬ 
ness  is  shown  in  Figure  28a.  Because  rutile  is  a  faster  material 
the  velocity  increases  with  layer  thickness  until  reaching  the 
slow  shear  velocity  of  the  substrate.  Surface  wave  coupling  is 
plotted  in  Figure  28b  for  two  transducer  configurations.  The  cou¬ 
pling  becomes  zero  for  layers  greater  than  10%  of  a  wavelength. 

Selected  numerical  solutions  for  rutile  layers  are  shown  in 
Figures  29-31.  The  displacements/potential  show  penetration  in¬ 
creasing  with  layer  thickness.  The  increased  penetration  with  in¬ 
creasing  layer  thickness  is  due  to  the  nearness  of  the  bulk  shear 
wave  in  the  substrate. 

Numerically  TiO  2/ST-Quartz  is  governed  by  equations  9  and  10 
of  section  3.2  which  contain  10  sub-solutions.  Displacement  and 
potential  vectors  associated  with  each  sub-solution  are  shown  to¬ 
gether  with  the  10  complex  weighting  factors  for  the  complete 
summation  solution  per  equation  5  of  Section  3.1. 
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4.1.6  Analysis  of  Al^ O^/ST-Quartz 

Layers  of  sapphire  (Al203)  on  ST-cut  Quartz  substrates  were 
analyzed.  Elastic  and  dielectric  constant  data  for  c-axis  orient¬ 
ed  sapphire  were  used  (Appendix  I) .  Surface  wave  velocity  vs  nor¬ 
malized  layer  thickness  is  shown  in  Figure  32a.  Because  sapphire 
is  a  fast  material  the  velocity  increases  with  layer  thickness 
until  reaching  the  slow  shear  velocity  of  the  substrate.  Surface 
wave  coupling  is  plotted  in  Figure  32b  for  two  transducer  config¬ 
urations.  The  coupling  becomes  zero  for  layers  greater  than  10% 

of  a  wavelength. 

Selected  numerical  solutions  for  sapphire  layers  are  shown 
in  Figures  33-35.  The  displacements  and  potential  show  penetrat¬ 
ion  increasing  with  layer  thickness.  The  increased  penetration 
with  increasing  layer  thickness  is  due  to  the  proximity  of  the 
bulk  shear  wave  in  the  substrate. 

Numerically  A^O^/ST-Quartz  is  governed  by  equations  9  and 
10  of  section  3.2  which  contain  10  sub-solutions.  Displacement 
and  potential  vectors  associated  with  each  sub-solution  are  shown 
together  with  the  10  complex  weighting  factors  for  the  complete 
summation  solution  per  equation  5  of  section  3.1. 
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Figure  34-  Layered  velocity  solution  for  A^O  3/ST-Quar  tz 

when  the  top  surface  is  shorted  and  the  interface 
is  unshorted.  Layer  thickness  is  10%  of  a 
wavelength . 
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Figure  35-  Layered  velocity  solution  for  Al£0 j/ST-Quartz 

when  the  top  and  interface  surface  are  shorted. 
Layer  thickness  is  10%  of  a  wavelength. 


4,2  Interfacial  Wave  Analyses 

In  this  section  highlights  of  a  study  of  Stoneley  waves  in 
anisotropic  crystals  are  presented.  For  illustration,  data  on 
copper  and  quartz  crystals  is  presented.  Copper  was  selected  to 
establish  a  basis  with  previous  publications.  Quartz  is  of  part¬ 
icular  interest  because  it  is  synonymous  with  frequency  control 
and  stability. 

The  approach  taken  in  this  interfacial  wave  study  is  shown 
in  Figure  36  where  a  view  along  the  Xj  direction  (normal  to  the 
X.*,  plane)  is  shown.  This  plane  is  called  the  interface  plane. 


Figure  36-  Coordinate  system  relative  to  the  interface  plane 

plane  for  Stoneley  wave  analysis.  Crystal  material 
parameters  are  rotated  with  respect  to  the  kx  dir¬ 
ection  along  which  the  interface  wave  propagates. 

In  order  to  obtain  interface  waves  in  a  single  material 
two  different  orientations  of  the  same  crystal  must  be  used.  As 
shown  in  Figure  36  each  orientation  (top  and  bottom  material)  is 
obtained  by  rotating  the  crystal  material  parameters  with  respect 
to  the  wave  propagation  vector.  Typically  the  top  material  is 
rotated  by  a  negative  angle  and  the  bottom  material  by  a  positive 
angle.  The  objective  of  our  study  was  to  determine  their  propert¬ 
ies  and  how  best  to  search  for  interface  waves. 
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4.2.1  Analysis  of  Copper 

The  interface  wave  velocity  as  a  function  of  material  rotat¬ 
ion  is  shown  in  Figure  37.  Unrotated  material  constants  for 


copper  (Appendix  I)  were  used.  The  top  and  bottom  material 


constants  were  obtained  by  rotating  negative  and  positive  amounts 


respectively. 
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Figure  37-  Stoneley  wave  velocity  vs  material  rotation  angl 
in  single  crystal  copper  (001  plane). 

Two  points  are  of  interest:  (1)  the  larger  angles  result  in  large 

bulk  like  shears  at  the  interface  and  (2)  angles  less  than  cut¬ 


off  do  not  possess  Stoneley  wave  solutions. 


A  Stoneley  wave  search  technique  was  developed  based  upon 

the  prediction  of  cut-off  using  techniques  of  Barnett  and  Lothe 
(Appendix  II) .  This  search  method  does  not  require  complete  sol¬ 
ution  of  the  wave  boundary  equations  and  greatly  simplified  the 


Stoneley  wave  searches. 


Stonely  wave  solutions  for  two  different  angles  are  shown  in 
Figure  38.  AT  i  10°  the  displacements  are  not  strongly  bound  to 
the  interface  however,  at  -  20°  the  wave  is  well  bound  and  decays 
within  3-5  wavelengths  of  the  interface. 

Stoneley  wave  solutions  were  found  to  be  strongly  influenced 
by  the  material  anisotropy  as  defined  by 

2  C44  /(C11  -°12  > 

Normal  copper  possesses  an  anisotropy  of  3.2  however,  artificial¬ 
ly  changing  C^to  obtain  an  anisotropy  of  1.52  and  the  effects 

this  has  on  Stoneley  wave  solutions  is  shown  in  Figure  39.  Lower 
anisotropy  lowers  the  cut-off  angle  and  causes  the  wave  to  be 
more  closely  bound  to  the  interface. 

4.2.2  Analysis  of  Quartz 

Using  the  above  search  techniques  (Appendix  II)  a  study  was 
performed  on  single  crystal  Quartz.  Although  these  waves  were  not 
found  to  exist  on  all  cuts  of  quartz  they  have  for  the  first  time 
been  found  in  quartz.  The  data  indicates  piezoelectrically  act¬ 
ive  and  well  bound  waves  are  possible.  This  result  suggests  a 
potentially  useful  resonator  structure  with  electrodes  encapsula¬ 
ted  at  the  material  interface. 

A  complete  listing  of  all  numerical  search  data  within  this 
report  is  not  possible  due  to  length  restrictions  hence  only  one 
Stoneley  wave  solution  for  quartz  will  be  presented.  Shown  in 

Figures  40  and  41  are  the  velocity  solutions  for  X-cut,  Z  propa¬ 
gating  Stoneley  waves  w?th  a  rotation  angle  of  -  20°.  The  effects 
of  shorting  the  electric  field  at  the  interface  are  shown  in  the 
detailed  plots  of  Figure  42. 


60 


RELATIVE  DISPLACEMENT/POTENTIAL  RELATIVE  DISPLACEMENT /POTENTIAL 


Figure  3 
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(b) 


8-  Stoneley  wave  solution  on  Copper  fjr  a 

rotation  angle  of  (a)  i  10°  and  (b)  -  20°. 


RELATIVE  DISPLRCEMENT/POTENTIRL  RELATIVE  DISPLACEMENT /POTENTIAL 


ROOT#  ]  +0.0000,-1.0225  Ul=  +1.000,  +0.000  BOT 

U2=  -.963,  +0.000  BOT 

U3”  +0.000,  -.014  BOT 

U4=+0.000E+00,+S.238E+ll  SUB 


CS 

cvi  w  ^  w  y 

W  *  *  *  * 

O  *  ^  i  H 

r  <-#  m  m  C 


?  *  *  * 
>  *  *  ¥ 

H 


c6  <r 


-  * -  * - 1 - 1 _ ± _ L _ i _ 

cncof^ujin^rcnM 


nuuN3xod/iN3waDundsia  3Aiitn3d 


CO  CO 

•  to  f-  H  3 

>  w  o  o  w 

>  CQ  CQ  I 

O  ( 

CO  CO  03  ( 

>  as  in  m  as 

*  +  in  ,h  + 

movinMi 
.  as  ess  •  •  as  f 

as  ®  <h  i  as  t 

as  •  +  as 

♦  as  •  < 

cs  +  as 

I  l 


i  to 

>  a-  a.  o  c 

»  o  o  &  < 

h  t*  i 
a+  to 

.  O  as  O 

i  t-.  as  ro  css  t-  c 

ca  to  +  c 

I  s®  H  U  O  t 

i  cs  •  «  as  as 

■  CS  ®  +  CM  ®  C 

>  •  +  ■ 


as  os  **r^  o  as  *as®a  »■  to  cs  cs  »cooa  *  as  m  as  *  as  in  as 
cancoooncots  cstocscscscntsocsr-Hcscs  cscoocscsascscs 
cavon  +  csio<,S4oiD's  +  a«)®  +  o«)®  +  ar^®  +  OMts  + 

CD*  •  CO  <5  •  •  CO  as  •  •  to  CS  «  •  to  cs  •  »  to  ©  •  •  to  CS  •  •  to 

•  I  nm  'inn  •  in  cs  css  •  +  o  a  •  t  ®  ®  *  cm  •  n  v  <h 

rH  I  M*  rH  +  V  H  M  +  a  H  +  a  H  +  O  H  I  I  M  H  I  +  CS 

+  cs  4-  cs  +  +  ®  +  cs  -f  cs  +  m  +  m 

CM  CM  O  CS  Cs  <-4  rH 

b  h  ii  +•  h  u  u  iMnn  +  niirt  +  Muii  +  iiuu  »  n  11  11  + 

»h<mcou  h(n  n  u  rH  cm  m  u  hoi  n  m  hnh  u  nwn  n  hcn  n  n 

tototoM-toOO^CSZJtovDtoto^rtototoM’atotoirtototoi*' 
to  to  to  to  to  to  3 


to  to  to  to 
as  oo  rM  cs 
Non® 
in  m  m  ® 

HHH® 
+  11  + 
*  *  *  * 
HHfO® 

cs  <s  cs  cs 
111  + 
to  to  to  to 
in  .h  in  3 
in  r-  in  ® 
C"  r-  t  <s 

NHHH 
111  + 
u  B  n  N 


+T  H  H  rl 

oooo 

I  I  I  I 

to  to  to  to 

+f  Q  l/l  V 

a>  as  in  in 
m  »-#  as 

m  as  m  as 
to  I  l  i  + 

££**•**■ 
o  v  O  H  N 
Hoaoo 
U  I  +  «  I 

<  to  to  to  to 

to  cm  r-  oo  co 

r~  M*  CS  rH 

U  n  ^  o  i/i 
z  •  •  •  • 

►-*  O'  rH  to  Os 

t-  I  +  I  + 

Z  0  I  II  H 

~  ~~  ~ 

mHMI/IP' 

u  ^  w  w  +-T 

2  to  to  to  to 


63 


Figure  AO-  Stoneley  wave  velocity  solution  for  XZ-Quartz 
when  the  interface  is  unshorted. 
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Figure  42-  Expanded  plot  of  the  displacements  and  potential 
of  Stoneley  waves  on  XZ-Quartz  when  the  interface 
is  (a)  unshorted  and  (b)  shorted. 


Stoneley  Interface-Wave  Accelerometer 
A  Stoneley  wave  is  a  type  of  elastic  wave  that  is  propagated 
along  the  interface  of.  two  elastic  solids.  Because  its  energy  re¬ 
mains  at  or  near  the  interface,  it  can  be  generated,  reflected 
and  detected  at  the  interface.  Since  Stoneley  wave,  derived  from 
surface  acoustic  wave  (SAW)  technology,  is  also  a  planar  techno¬ 
logy,  and  hence  can  be  fabricated  by  photolithographic  techni¬ 
ques.  That  is  why  Stoneley  wave  devices  can  be  fabricated  at 
costs  comparable  to  those  for  integrated  circuits. 

The  interface  Stoneley  wave  resonator  consists  of  a  pair  of 
reflective  gratings  on  a  surface,  defining  a  resonant  cavity  for 
interface  acoustic  waves.  These  waves  can  be  generated  and 
detected  by  an  interdigital  transducer,  which  is  piezoelectr ical- 
ly  coupled  to  the  interface  strain  variations  of  the  acoustic 
waves.  This  resonant  cavity,  shown  in  Fig. 43,  is  coupled  to  a 
feedback  amplifier  to  maintain  the  energy  in  the  cavity. 


Figure  43-  Stoneley  Interface  Wave  Configuration 
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The  bit  quantization  for  an  interface-wave  accelerometer  is 
in  the  same  range  as  those  of  high  precision  accelerometers.  For 
an  interface  wave  accelerometer  with  320  MHz  frequency  and  proof 


mass  loading  to  produce  200  ppm  strain  at  20g  input  acceleration, 
the  bit  quantization  will  be: 


(20  g)  (32  fps/sec/g) 


6  **6 
(320  x  10  bit/sec)  (200  x  10  p/p) 


=  0.01  fps/bit 


The  short-term  stability  of  the  interface  wave  resonator  (as 


suming  for  the  worst  case  is  similar  to  surface  acoustic  wave)  is 


in  the  order  of  0.1  parts  per  billion. 


The  accelerometer  overall  design  has  a  direct  bearing  on  the 
complexity  of  the  inertial  guidance  system.  The  digital  nature  of 
the  interface  wave  accelerometer  output  results  in  an  extremely 
simple  system  architecture.  Shown  in  Fig. 44  is  a  baseline  guidan¬ 
ce  system  block  diagram  which  shows  how  the  interface  wave  accel¬ 
erometer  signals  are  processed  to  generate  not  only  acceleration 
but  also  velocity  and  displacement  signals  (words)  for  the  gui¬ 
dance  system. 
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OSCILLATOR  No.  1 


Figure  44-  Proposed  Baseline  Interface  Wave  Accelerometer  in 

Inertial  Guidance  System 


The  mixer  output  is  a  frequency  proportional  to  acceleration 
and  the  acceleration  counter  generates  a  binary  count  with  each 
clock  reset  from  the  microprocessor.  The  output  of  the  accelerat¬ 
ion  counter  is  summed  (integrated)  in  another  counter  which  then 
generates  a  binary  word  representing  velocity.  Displacement  is 
computed  by  summing  the  velocity  counter  output  and  transmitting 
the  displacement  word  to  the  microprocessor .  The  microprocessor 
provides  the  appropriate  scaling  and  generates  the  needed  guidan¬ 
ce  and  control  signals. 

,  The  digital  system  proposed  is  capable  of  providing  acceler¬ 
ation  signals  with  more  than  12  bits  accuracy  due  to  the  large 
dynamic  range  of  such  proposed  accelerometer .  The  system  is 
unique  in  that  12  bit  accuracy  is  provided  without  using  analog- 


to-digitai  converters.  This,  in  turn*  will  enable  the  design  of  a 
system  with  minimum  power  requirements.  Starting  with  this 
baseline  system#  the  detailed  design  will  be  the  subject  for  the 
proposed  study. 
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4.4  Comparison  With  Experimental  Results 

The  theoretical  results  of  Section  4.1  were  compared  with 
experimentally  published  results.  The  experimental  method  con¬ 
sisted  of  monitoring  the  frequency  (proportional  to  velocity) , 
resonant  resistance,  and  quality  factor  (Q)  of  SAW  resonators 
while  a  crystalline  film  was  deposited  (sputtered)  under  vacuum. 
The  following  films  were  evaluated  by  this  method:  Mgo,  SiO  , 

X 

Y  ^0  y  and  AIN.  The  experimental  results  are  shown  in  Figures 
45  and  46. 

Because  the  experimental  results  are  based  upon  time  of  dep¬ 
osition  and  not  absolute  film  thickness  it  was  difficult  to  com¬ 
pare  directly  theory  and  experiment.  Nevertheless  it  can  be  ass¬ 
umed  that  thickness  is  proportional  to  time  for  discussion  purp¬ 
oses.  Comparing  the  theoretical  results  (MgO)  of  Figure  9  with 
Figure  45,  several  points  are  noted:  (1)  the  velocity  or  frequency 
increases  with  film  thickness  and  (2)  after  an  initial  drop  in  Q 
the  loss  is  relatively  constant  with  only  a  small  increase  with 
thickness.  The  increase  in  frequency  with  thickness  is  predicted 
because  we  believe  the  film  grown  was  considerably  less  than  10% 
of  a  wavelength. 

Experimental  results  for  SiOx  films  are  compared  to  theore¬ 
tical  predictions  using  Figure  18.  The  important  points  are:  (1) 
experimental  frequency  drops  with  film  thickness  in  contrast  to 
theory  which  predicts  an  increase  with  film  thickness  and  (2)  ex¬ 
perimental  results  indicate  a  sharp  increase  in  loss  at  an  inter¬ 
mediate  thickness.  The  later  is  compatible  with  what  is  predict- 
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ed  for  film  thickness  in  the  range  than  5-10%  of  a  wavelength. 
The  discrepency  between  the  frequency  behavior  is  explained  as 
follows:  For  film  thickness  less  than  5%  of  a  wavelength  the  fre¬ 
quency  changes  by  a  very  small  amount,  however  this  is  masked  by 
a  relatively  large  temperature  increase  during  sputtering.  Half 
way  into  the  film  deposition  there  is  a  decrease  in  the  rate  of 
frequency  drop  and  at  the  same  time  the  loss  increases  dramatic¬ 
ally.  We  believe  this  is  where  the  film  thickness  is  approaching 
5-10%  of  a  wavelength  resulting  in  increased  wave  penetration  and 
a  sharp  drop  in  coupling. 

Experimentally  films  of  O3  were  the  most  well  behaved  with 
a  uniform  drop  in  frequency  during  deposition,  a  relatively  small 
increase  in  resistance,  and  a  drop  in  Q.  The  decrease  in  freque¬ 
ncy  and  the  low  loss  was  predicted  because  the  surface  wave  does 
not  couple  into  a  bulk  wave  for  Y^ 03  layers  at  any  thickness. 

Exper imentaly  films  of  AIN  showed  an  increase  in  frequency 
and  a  uniform  increase  in  loss  during  deposition.  Near  the  end 
of  the  deposition  the  resistance  increase  was  accelerated.  Theor¬ 
etically  the  increase  in  frequency  was  predicted.  The  large  inc¬ 
rease  in  resonator  resistance  and  loss  near  the  end  of  the  depos- 
iton  indicates  the  film  was  near  5-10%  of  a  wavelength  hence  the 
piezoelectric  coupling  drops  and  coupling  to  a  bulk  wave  occurs. 


72 


5.0  CONCLUSIONS  AND  RECOMMENDATIONS 


This  study  has  established  the  theoretical  foundations  for 
using  higher  order  Rayleigh  and  interfacial  (Stoneley)  waves  for 
acoustic  devices  such  as  SAW  resonators  and  sensors.  Techniques 
for  analyzing  layered  SAW  devices  as  well  as  multi-layered  inter¬ 
facial  waves  were  developed  and  applied  to  many  anisotropic  mat¬ 
erials  as  part  of  this  study. 

Layered  SAW  device  performance  using  materials  such  as  MgO, 
SiOx,  AIN,  ¥203,  A1203,  and  Ti02  was  theoretically  analyzed  and 
the  results  compared  with  experimental  results.  Shifts  in  freque¬ 
ncy,  loss  to  bulk  mode  coupling,  and  changes  in  piezoelectric 
coupling  can  now  be  predicted  for  any  layer  material.  From  this 
study  it  is  concluded  that  passivation  ,  for  the  first  order 
Rayleigh  mode  cannot  be  achieved  using  layers  of  material.  This 
is  due  to  the  fact  that  the  particle  displacements  of  the  surface 
wave  are  at  the  surface  of  any  deposited  layer  and  hence  a  decre¬ 
ase  in  piezoelectric  coupling  because  of  the  separation  between 
the  elastic  displacements  of  the  wave  and  the  potential  generat¬ 
ing  electrodes  at  the  interface. 

For  future  studies  of  layered  SAW  devices  it  is  recommended 
that  temperature  effects  and  power  flow  calculations  be  added  to 
the  existing  analysis  techniques.  These  properties  can  be  used 
to  predict  temperature  compensated  materials  and  determine  orien¬ 
tations  where  power  flow  is  co-linear  with  the  wave  vector. 
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An  analysis  of  interfacial  waves  in  multi-layered  media  was 
completed  during  this  phase  I  study.  A  method  of  analyzing  inter¬ 
facial  waves  was  developed  which  included  search  methods  which 
greatly  reduced  the  time  needed  to  determine  the  existence  of 
Stoneley  waves.  Using  these  methods,  Stoneley  waves  were  for  the 
first  time  predicted  in  single  crystal  quartz.  Stoneley  waves 
which  are  well  bound  and  piezoelectr ically  active  were  found  to 
exist  for  several  orientations  rf  quartz.  This  is  an  important 
discovery  since  it  may  result  in  a  new  class  of  acoustic  resonat¬ 
ors  and  sensors  with  stability  superior  to  existing  SAW  devices. 
Unlike  layered  SAW  devices,  Stoneley  wave  structures  are  well 
passivated  and  require  no  hermetic  enclosure;  thus  they  will  have 
better  stability  and  lower  cost. 

Based  upon  the  results  of  this  phase  I  study,  recommendat¬ 
ions  for  future  work  are: 

1.  To  add  temperature  sensitivity  and  power  flow  analysis 
to  the  existing  theoeretical  techniques. 

2.  To  design,  fabricate,  and  test  interface  wave  resonators 
and  to  compare  their  characteristics  with  theoretical 
predictions . 

3.  To  design,  fabricate,  and  test  interface  wave  sensors  as 
inertial  navigation  (acceleration)  sensors. 
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7.0  APPENDICES 


THE  MATERIAL  IS  ST-CUT  QUARTZ 

+8.670E+10  +2.749E+10  -8.587E+09 
+2.749E+10  +9.6S3E+X0  -4.830E+09 
-8. 587E+09  -4.830E+09  +1.306E+11 
-1.038E+09  -1.339E+10  +1.872E+0! 
+0 . 000E+00  +0.000E+00  +0 . 000B+00 
+0.000E+00  +0.000E+00  +0.000E+00 

+1.710E-01  -3 . 832E-02  -1.327E-01 
+0 . 000E+00  +0 . 000E+00  +0.000E+00 
+0.000E+00  +0 . 000E+00  +0.000E+00 

+3 . 920E-1 1  +0 . 000E+O0  +0 . 0G0E+00 
+0.000E+00  +4.017  E-ll  -8.972E-13 
+0 . 000E+00  -8 . 972E-13  +4.003B-11 

DENSITY  =  2650 


■iriE  MATERIAL  IS  FUSED  QUARTZ 

+7 . 850E' ' 0  +1.6.3E+10  +1.610E+10 
+1.610E*-0  +7 . 850E+10  +1.610E+10 
+1 . 610E+10  +1.610E+10  +7.850E+10 
+0 . 000E+00  +0.000E+00  +0 . 000E+00 
+0 . 000E+00  +0.000E+00  +0.000E+00 
+0 . 000E+00  +0.000E+0O  +0 . 000E+00 

+0. 000E+00  +0.000E+00  +0.000E+00 
+0 . 000E+00  +0.000E+C0  +0.000E+00 
+0 . 000E+00  +0.000E+00  +0 . 000E+00 

+  3 . 320E-11  +0.000E+00  +0.000F.+00 
+0.000E+00  +3.320E-11  +0.000E+00 
+0.WJ0E+00  +0 . 000E+00  +3 . 320  E-ll 

DENSITY  =  2200 


-1.038E+09  +0.000E+00  +0 . 000E+00 
-1 . 339E+10  +0.000E+00  +0.0O0E+00 
+1 , 872E+09  +0.000E+00  +0.000E+00 
+  4 , 1 17  E+10  +0.000E+00  +0.000E+00 
+0 . 000E+00  +3 . 032E+10  -7 . 593E+09 
+0 . 000E+00  -7.593E+09  +6.743E+10 

-8.205E-02  +0.000E+00  +0.000E+00 
+  0 . 000  £+00  -6 . 653E-02  -9.903E-02 
+0.000E+00  -7 . 197E-02  -1 . 071E-01 


+0 . 000E+00  +0.000E+00  +0.000E+00 

+0.000E+00  +0.000E+00  +0.000E+00 

+0.000E+00  H0.000E+00  +0.000E+00 

+  3 . 120  E  +  10  +0.000E+00  +0.000E+00 

+  0.000E+00  +3 . 120E  +  10  +0 . 000E+00 

+0 . 000E+00  +0.000EV00  +3 . 120E+10 

+0.000B+00  +0.000E+00  +0.000E+00 

+0.000E+00  +0.000E+00  +0.000E+00 

+0 . 000E+00  +0.000E+00  +0.000E+00 
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THE  MATERIAL  IS  YITTRIUM  OXIDE 


+  2 . 270E+11  +I.  380E+11  +  1.380B-I11 

+1.380E+11  +2.270E+11  +1.380E+11 

+1.380E+11  +1.380E+11  +2 . 270E+11 

+0 . 000E+00  +0.000E+00  +0 . 000E+00 

+0.000E+00  +0.000E+00  +0.000E+20 

+0.000 E+00  +0.000E+00  +0.000E+00 

+0 . 000E+00  +0.000E+00  +0 . 000  E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000  E+00 

+0 . 000E+00  +0 . 000E+00  +0 . 000  E+00 

+  1 . 151E-10  +0 . 000  E+00  +0 . 000  E+00 

+0 . 000E+00  +1 . 15IE-10  +0 . 000  E+00 

+0. 000E+00  +0.000E+00  +1 . I51E-10 

DENSITY  =  5010 


THE  MATERIAL  IS  UNROTATED  COPPER 

+  1 • 710E+1 1  +1 . 239E+11  +1 . 239  E+ll 

+  1 . 239E+11  +1 . 710E+11  +1.239  E+ll 

+1 . 239E+11  +1 . 239E+11  +1.719E+11 

+0 . 000  E+00  +0.000E+00  +0 . 000  E  +  00 

+0 . 000E+00  +0 . 000E+00  +0 . 000  E+00 

+  0 . 000  E+00  +0 . 000  E+00  +0 . 000E+00 

+0 . 000E+00  +0 . 000  E»00  +0.000E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000E+00 

+0.000E+00  +0 . 000  E+00  +0 . 000E+00 

+0.000E+00  +0 . 000E+00  +0 . 000E+00 

+0.000E+00  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000E+00  +0 . 000  E+00  +0.000E+00 

DENSITY  =  8950 


THE  MATERIAL  IS  RUTILE  (031) 

+2 . 660E+11  +1.730E+11  +1.360E+11 

+1.730E+11  +2 . 663  E+ll  +1.360E+11 

+1 , 360E+11  +1.360E+11  +4.700E+11 

+0.0003+00  +0 . 000E+00  +0 . 000E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000E+00  +0 . 000  E+00  +0 . 000E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000  E  t-00 

+0 . 000E+00  +0 . 000E+00  +0 . 000  E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000E+&0 

+7.610E-10  +0 . 000  E+00  +0 . 000E+00 

in  n  n  /i  *■»  •  n  n  .  ri  n  •  l  ft  ■  si  An/n<  ,  a,i 
<tf«t/UUL  TVV  T  /  «  OlDL“ll/  lU.WUl/iilUl/ 

+0 . 000E+00  +0.000E+00  +1 . 504E-09 

DENSITY  =  4280 


+0 . 000  E+00  +0 . 000E+00  +0 . 000E+00 

+  0 . 003  E+00  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000  E+00  +0 . 000  E+00  +0 . 000E+00 

+  6 . 860E+10  +0.000E+00  +0 . 000  E+00 

+  0 . 000  E+00  +6.860E+10  +0 . 000E+00 

+  0 . 000  E+00  +0 . 000  E+00  +6 . 860E+10 

+0 . 000  E+00  +0.000E+00  +0 . 000  E+00 

+  0 . 000  E+00  +0 . 000  E+00  +0 . 000  E+00 

+0 . 000  E+00  +0 . 000  E+00  +0.000E+00 


+  0 . 000  E+00  +0.000E+00  +0 . 000E+00 

+  0 . 000E+00  +0 . 000  E+00  +0 . 000  E+00 

+0 . 000E+00  +0 . 000  E+00  +0 . 000  E+00 

+  7 . 5G0E+10  +0 . 000  E+00  +0 . 000  E+00 

+0 . 000  E+00  +7.560E+10  *0 . 000  E+00 

+0 . 000E+00  +0 .000  E+00  +7 . 560E+10 

+0 . 000E+00  +0 .000  E+00  +0 . 000E+00 

+0.000E+00  +0.000E+00  +0 . 000  E+00 

+  0 . 000  E+00  +0 . 000E+00  +0 . 000  E+00 


+  0.0UIIE+00  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000E+00  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000  E+00  +0 . 000  E+00  +0 . 000E+00 

+  1 . 240  E+ll  +0 . 000  E+00  +0 . 000  E+00 

+  0 . 000  E+00  +1 . 24  0  E+ll  +0 . 000  E+00 

+0 . 000 E+00  +0 . 000 E+00  +1 . 890E+11 

+  0 . 000  E+00  +0 . 000  E+00  +0 . 000  E+00 

+  0.000E+00  +0 . 000  E+00  +0 . 000E+00 

+  0 . 000E+00  +0 . 000  E+00  +0 . 000E+00 
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THE  MATERIAL  IS  Z_CUT_X_PROP_MgO 


+? . 964E+1I 
+9.G00E+10 
+9.500E+10 
+0.000E+00 
+0 . 000E+00 
+  0 . 000E+00 

+0 . 000  E+00 
+0 . 000E+00 
+  0 . 000  E+00 

+8.544E-11 
+0.000 E+00 
+  0 . 000E+00 


+9. 500E+10 
+2.964E+1I 
+9.500E+10 
+0 . 000E+00 
+0 . 000  E+00 
+0 . 000  E+00 

+0 . 000E+00 
+  0 . 000E+00 
+  0 . 000E+00 

+0 . 000  E+00 
+8.544E-1I 
+0.000E+00 


+9 . 500E+10 
+9 . S00E+10 
+2.964E+11 
+  0 . 000  E+00 
+  0 . 000  E+00 
+  0 . 000  E+00 

+0 . 000  E+00 
+0 .000E+00 
+0 . 000  E+00 

+0 . 000E+00 
+0 . 000E+00 
+8.544E-11 


+0.H0HB+00 
+  0 . 000  E+00 
+0.000E+00 
+  U560E+1I 
+0 . 000E+00 
+  0 . 000  E+00 

+  0 . 000  E+00 
+0.000E+00 
+  0 . 000  E+00 


+0,000 E+00 
■i  0  . 000  E+00 
+0 . 000E+00 
+  0 .000  E+00 
+1 . 560E+11 
+  0 . 000  E+00 

+0 . 000  E+00 
+0 . 000  E+00 
+  0 . 000  E+00 


DENSITY  =  3583 


THE  MATERIAL  IS  SAPPHIRE  (001) 


+4 . 970E+11 
+1 . 640E+11 
+1. 110E+11 
-2 . 350E+10 
+0 . 000  E+00 
+0 . 000E+00 

+0 . 000E+08 
+0 . 0O0E+00 
+0 . 000E+00 

+  8 . 280E-11 
+0 . 000E+00 
+0 . 000E+00 


+1 . 640E+11 
+4.970E+11 
+1 . 110E+11 
+2.350E+10 
+  0 . 000  E+00 
+0 . 000  E+00 

+0.000E+00 
+  0 . 000  E+00 
+0.000E+00 

+0.000E+00 
+  8.280E-U 
+  0 . 000 E+00 


+  1 .110E+11 
+  1.U0E+11 
+4 .980E+11 
+0 . 000E+00 
+0 . 000  E+00 
+0 . 000  E+00 

+  0 . 000  E+00 
+  0 . 000  E+00 
+0 . 000E+00 

+0 . 000E+00 
+  0 . 000E+00 
+  1 . 020E-10 


-2 . 350E+10 
+2 . 350E+10 
+0 . 000  E+00 
+  1 . 470E+11 
+  0 . 000  E+00 
+  0 . 000  E+00 

+0.000 E+00 
+0 . 000E+00 
+0 . 000  E+00 


+0.000E+00 
+0.000E+00 
t0 . 000  E+00 
+  0 . 000  E+00 
+  1 . 470E+11 
-2 . 350E+10 

+0 . 000E+00 
+0.000E+00 
+  0 . 000  E+00 


DENSITY  =  3980 


THE  MATERIAL  IS  ZX  AL  NITRIDE 


+  1.2O0E+U 
+1 . 20OE+11 
+3 . 950E+1 1 
+0 . 000E+00 
+0 . 000E+00 
+  0 . 000  E+00 


+0.000E+00 
+0.000 E+00 
+  0 . 000  E+00 
+1 . 180E+11 
+  0 . 000  E+00 
+0 . 000E+00 


+0 . 000E+00 
+0 . 000E+00 
+0 .000E+00 
+  0 . 000E+00 
+  1 . 180E+11 
+0 . 000E+00 


+3 . 450E+11 
+1.250E+11 
+1 . 200E+11 
+0 . 000  E+00 
+0 . 000E+00 
+0 . 000E+00 

+0 • 000E+00 
+0 . 000E+00 
-5. 800E-01 

+  8 . 000E-11 
+0 . 000  E+00 
+  0  •  000  E+00 


+1.250E+1L 
+3 . 450E+1 1 
+  1 . 200E  +  11 
+  0 .000E+00 
+  0 . 000E+00 
+  0 . 000E+00 

+0 . 000E+00 
+  0 . 000E+00 
-5 . 800E-01 

+0.000E+00 
+8.000E-11 
+0 . 000E+03 


+0 . 000E+00 
+0.000E+00 
+1 . 550E+00 

+0 . 00UE+00 
+0 . 000E+00 

-  r. 

t  y  *  Ovv  r 


+0 . 000E+00 
-4 . 800E-0 l 
+0 . 000  E+00 


-4 . 800E-01 
+0 . 000E+00 
+0 . 000E+00 


DENSITY  =  3260 


+0 . 000E+00 
+  0 , 000  E+00 
+0 .000E+00 
+0 .000E+00 
+0 . 000S+00 
+1 . 560E+11 

+  0 . 000  E+00 
+  0 . 000  E+00 
+0 . 000  E+00 


+0.000E+00 
+  0 . 000  E+00 
+0 . 000E+00 
+  0 . 000  E+00 
-2 . 350E+10 
+1 . 665E+11 

+0 . 000  E+00 
+0.000E+00 
+0.000E+00 


+0 , 000E+00 
+0.00OE+00 
+  0 . 000  E+00 
+0 . 000E+00 
+  0 . 000  E+00 
+1 . 100E+11 

+  0 . 000E+00 
+0.000E+00 
+  0 . 000E+00 
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7.2  Appendix  II 


The  basis  of  a  search  technique  for  Stoneley 
ained  in  the  paper  reprinted  in  this  appendix, 
of  formulating  the  wave  impedance  is  discussed 
detail  by  the  authors.  Using  these  techniques  a 
for  Stoneley  waves  in  piezoelectric  materials  was 


waves  is  cont- 
The  importance 
in  considerable 
search  method 
developed . 
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On  the  existence  of  surface-wave  solutions  for  anisotropic 
elastic  half-spaces  with  free  surface 

J.  Lothe 

Institute  of  Physics.  University  of  Oslo,  Oslo,  Norway 

D.  M.  Barnett 

Department  of  Materials  Science  and  Engineering,  Stanford  University,  Stanford,  California  94305 


A  proof  is  developed  that  for  a  given  direction  of  propagation  on  the  free  surface  of  a  half-infinite 
anisotropic  crystal,  a  surface-wave  solution  with  a  certain  phase  velocity  v*  <  vt,  where  vt  is  the  limiting 
velocity,  will  always  exist,  except  in  the  special  case  when  the  bulk  wave  defining  the  limiting  velocity 
satisfies  the  condition  of  a  free  surface.  The  proof  is  in  terms  of  the  surface  impedance,  which  relates  the 
amplitude  at  the  surface  of  a  surface  wave  with  the  external  forces  needed  at  the  surface.  The  properties  of 
the  impedance  as  a  function  of  phase  velocity  determines  whether  a  surface  wave  not  requiring  external 
forces  at  the  surface  exists  for  a  certain  phase  velocity.  The  proof  is  valid  also  in  the  case  of  degeneracies 
in  the  eigenvalue  problem  entering  the  analysis. 


I.  INTRODUCTION 

Recently  the  present  authors'  provided  a  proof  for  the 
existence  of  surface-wave  solutions  for  elastic  half¬ 
spaces  with  a  free  boundary.  However,  the  proof  was  in 
terms  of  the  Lagrangian  for  moving  dislocations,  a 
quantity  not  well  known  among  surface- wave  re¬ 
searchers,  Here  an  alternative  proof  in  terms  of  the  so- 
called  surface-wave  impedance  is  given.  Usually  a  sur¬ 
face  wave  with  given  phase  velocity  and  given  amplitude 
at  the  surface  requires  external  forces  at  the  surface. 
The  surface  impedance  relates  amplitude  at  the  surface 
with  external  forces  needed  at  the  surface.  However, 
when  the  impedance  satisfies  certain  requirements, 
surface  waves  which  do  not  require  external  forces  at 
the  surface  may  exist.  In  this  paper  we  will  analyze  the 
impedance  as  a  function  of  velocity  to  see  whether  for 
certain  phase  velocities  the  impedance  satisfies  the  re¬ 
quirements  for  surface  waves  in  the  case  of  free  sur¬ 
faces. 

In  Secs.  II  and  III  an  integral  form  of  the  general 
theory  for  surface  waves  maintained  by  external  forces 
is  developed  and,  in  particular,  a  convenient  integral 
expression  for  the  impedance  is  derived.  As  in  our 
previous  paper*  the  analysis  presupposes  no  degeneracy 
in  the  eigenvalue  problem  appearing.  In  Sec.  IV  the 
complications  that  arise  with  degeneracy  are  considered 
and  it  is  concluded  that  the  integral  expression  for  the 
impedance  remains  valid.  In  Sec.  V  the  properties  of 
the  impedance  are  listed.  Finally,  in  Sec.  VI  we  turn  to 
the  special  problem  of  surface- wave  propagation  on  free 
surfaces.  This  is  the  real  problem  to  which  we  are 
addressing  ourselves  in  this  paper.  It  will  be  shown  how 
this  problem  can  be  discussed  in  terms  of  the  properties 
of-the  impedance  developed  in  Secs.  Ill— V.  The  first 
sections  of  the  paper  are  necessary  prerequisites  for 
the  discussion  in  Sec.  VI. 

The  exposition  in  what  follows  will  be  kept  quite  brief 
by  appropriate  reference  to  earlier  work. 

II.  INTRODUCTORY  THEORY 

The  surface  has  normal  n°,  see  Fig.  1.  We  consider 


harmonic  waves  traveling  along  the  surface  in  the  m° 
direction 

u(x)  -  A(n°x)  expf/A’(m"x-  !•/)],  (2. 1) 

where  r  is  the  phase  velocity  for  propagation  along  the 
surface,  x  is  the  position  vector  with  a  point  on  the  sur¬ 
face  as  origin,  and  n°x  is  thus  the  depth  below  the  sur¬ 
face.  With  respect  to  a  fixed  external  Cartesian  lab¬ 
oratory  set  of  unit  vectors,  x  has  components  x,,  x2, 
and  x3. 

The  amplitude  A  as  a  function  of  depth  depends  on  the 
boundary  conditions  at  the  surface  and  the  boundary  con¬ 
ditions  at  infinite  depth.  For  given  k,  m,n,  and  v,  sev¬ 
eral  independent  waves  of  the  type  (2. 1)  will  be  possible. 
These  are  derived  in  Secs.  Ill  and  IV. 

Referring  also  the  elastic  coefficients  to  the  external 
laboratory  system,  the  elastic  force  per  unit  area 
parallel  with  the  surface  is 

■—  (2.2) 

or,  by  Eq.  (2. 1) 

Jj  —  il;()i<lCIJk,»ini)Ab  +  i  (2.3) 

where 

y  =  n°x  (2.4) 

is  a  convenient  notation  for  the  depth  below  the  surface. 
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FIG.  1.  Coordinate  systems  for  the  description  of  surface 
waves  In  a  half-infinite  crystal,  n°  is  the  ti.ner  unit  normal. 


If  an  observer  travels  with  a  wave  with  the  phase 
velocity  v  of  the  wave,  he  sees  a  stationary  wave  field, 
and  the  condition  for  dynamic  equilibrium  can  be  stated 
similarly  to  the  condition  for  equilibrium  in  the  static 
case  as 

c"“1^7=0'  <2-5> 

where  CUkl  are  the  dynamic  elastic  coefficients 

CiJkl  —Cu*i  ~  (2. 6) 

Introducing  the  vectors 

L  =  »rlJ,  (2.7) 

the  equilibrium  condition  (2. 5)  can  be  expressed  in  the 
form 

(2.8) 

where  5  is  the  six-dimensional  vector 

{=(£)  (2.9) 

and  N4  is  the  6X6  matrix 

/  (w°«)*1  («°m°)  (u°m°)-1  \ 

N0=- 

\(MiOnO)(jiO(!O)‘10iO»i°)  -  (»1°»H0)  (wi°«°)(»*0«0)“1/ 

(2. 10) 

The  3x3  matrices  ( ab )  occurring-  in  the  block  rep¬ 
resentation  of  N°  are  defined  by 

(ab)jk=alC{jkibi.  (2,11) 

The  symmetry  of  C'IJk,  is  such  that 
(1 ab)  =  (ba)T , 

where  T  means  transposition. 

The  above  is  essentially  the  formalism  of  Stroh2  and 
Ingebrigtsen  and  Tonning. 3 

In  terms  of  the  above  notation,  the  relation  between 
L  and  A  can  be  written 

L  =  -  (»iV)A+jfc-'(«V)-^- A  .  (2.12) 

Now  formally  construct  the  operator  N  by  the  pre¬ 
scription  (2. 10)  in  terms  of  a  unit  vector  system  (m,n) 
rotated  by  <p  relative  to  the  system  (m0,n°).  However, 
note  that  m°  in  the  definition  (2. 6)  of  C'Uk,  is  kept  fixed 
as  parallel  with  the  surface  and  the  direction  of  wave 
propagation.  The  resulting  N  is  a  function  of  <p,  N(p), 
and  as  shown  by  Barnett  and  Lothe4 


i.e. ,  ^ith  v  =  0  in  Eq.  (2.6),  (mm)  is  positive  definite 
by  elastic  stability, 

ll(m»i)||^0,  '•  v  =  0.  (2.16) 

Thus 

||mw||>0,  *» < >  all  0  (2.17) 

and  likewise 

||(mi)!l>0,  all  <p  when  v<vL 
so  that  (««)"'  and  also  N(<p)  exist  for  all  <p  when 

The  preceding  is  general.  Now  we  will  discuss  sep¬ 
arately  the  case  of  no  degeneracy  and  degeneracy  and 
show  that  the  same  integral  expression  for  the  surface- 
wave  impedance  applies  in  either  case. 

III.  NO  DEGENERACY 

Assume  that  A(n°x)=A(y)  in  Eq.  (2. 1)  is  of  the  form 

A(y)  =  A„  exp (ikp^y).  (3.1) 

Then  Eq.  (2. 8)  leads  to  the  eigenvalue  problem 

N“  «„=/>“«„,  (3.2) 

which  shows  that  the  p°a  must  be  the  roots  of 

IIN°  - />al  11  =  0.  (3.3) 

This  is  a  sixth-order  equation  with  real  coefficients, 
so  there  will  be  six  roots  occurring  in  pairs  of  complex 
conjugates.  Real  roots  cannot  occur  for  v  <  vL.  \\.v=vL, 
one  pair  of  roots  coalesces  into  one  real  root. 

When  all  pa  are  different  (no  degeneracy),  it  can  be 
asserted  by  general  theorems  (Nering5  and  Pease4)  that 

the  six  form  a  complete  set,  (3.4) 

or  said  in  a  more  detailed  fashion,  the  six  5a  are  linear¬ 
ly  independent  and  any  vector  5  in  the  six-dimensional 
vector  space  can  be  represented  by  a  linear  combination 
of  the  six  50 .  This  explains  the  meaning  of  the  word 
"complete”  in  theorem  (3.4). 

Tims  in  the  case  of  all  pa  different,  the  general  solu¬ 
tion  of  Eq.  (2. 8)  can  indeed  be  constructed  by  a  super¬ 
position  of  partial  solutions  of  the  type  (3. 1). 

Now  consider  the  eigenvalue  problem 

N(</>)^a=/,a(,/')^a,  (3.5) 

From  Eq.  (2. 13)  it  follows  that  the  eigenvectors  for 
Eq.  (3.5)  are  the  same  as  the  5a  of  Eq.  (3;  2).  The 
eigenvectors  5a  are  invariant  to  rotation. 

The  eigenvalues  must  vary  with  <p  as 


^■N  =  -(I+N2).  (2.13) 

The  operator  N(0)  exists  for  all  <p  up  to  the  limiting 
velocity  Dj,.  The  limiting  velocity  is  the  lowest  velocity 
for  which  there  exists  an  orientation  for  the  unit  vector 
m,  rax,  such  that 

II(wi»»l)II  =  0,  (2.14) 

allowing  a  bulk  wave  of  amplitude  A L  satisfying 

(mx»jx)Ax  =  0  (2.15) 

and  traveling  in  the  direction.  In  the  static  case, 


+Pl)  (3.0 

according  to  Eq.  (2.13).  From  Eq.  (3.6)  it  follows  that 


+  or  -  according  to  whether  we  consider  roots  with  a 
positive  or  negative  imaginary  part. 


Introducing 


(3.8) 
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where 


S=--~  f  (mm)'1  (nm)d<p, 

•'q 

(3.9) 

Q=-  -A-  f  (« n)’1  dip, 

■'0 

(3. 10) 

B=-  j  [(»»»)(««)*'  (/mi)- (mm)] d<p, 

(3.11) 

0 

it  now  follows  that 

Hi. 

(3. 12) 

The  first  three  of  these  equations  are 

SA„  +  QL0  =  ±  i'Aa . 

(3. 13) 

For  v<vL,  (mm)'1  is  positive  definite  and  consequently 
Q  is  negative  definite.  Thus  Q"1  exists  and 

L«=Q*‘(±»I-S)A0.  (3.14) 

Let  a  =  1,2, 3  be  the  solutions  with  positive  imaginary 
part  in  p„.  These  are  the  solutions  which  show  expo¬ 
nential  decay  with  increasing  y  and  which  are  to  be  used 
in  the  construction  of  surface  waves.  For  these  solutions 

L«=Q"‘(*I-S)Aoi  a  =  1,2, 3.  (3.15) 

Now  suppose  that  the  three  Aa  are  linearly  dependent. 
But  then  by  Eq.  (3. 15),  the  same  linear  dependence 
would  exist  between  the  corresponding  L„  and  conse¬ 
quently  between  the  corresponding  in  contradiction 
with  theorem  (3.4). 

Thus  the  three  Aa,  a  =  1,2, 3,  form  a  complete  set  in 
three-dimensional  space  for  all  v  <  vL.  The  same  con¬ 
clusion  was  reached  by  Stroh7  by  more  intricate 
reasoning. 

Since  A„,  £*  =  1,2,3,  form  a  complete  set,  an  ar¬ 
bitrary  amplitude  A  on  the  surface  can  be  obtained  by 
superposition  of  exponentially  decaying  solutions  pro¬ 
vided  a  force  amplitude 

J  =  -  i/fLi = At(Q_1  +  jQ’*S)A  (3.1G) 

is  maintained  at  the  surface.  We  can  write 

J--ZA,  (3.17) 

where 

Z=-A(Q-‘  +iQ"'S)  (3.18) 


where  I  is  block  notation  for  a  three-dimensional  unit 


matrix. 

Written  out  in  components  Eq.  (3. 21)  yields 

QB+SJ  +  1  =  0,  (3.22) 

SQ+QSr  =  0,  (3.23) 

BS  +  SrB  =  0,  (3.24) 

BQ  +  (Sr)2  +  I  =  0.  (3.25) 

From  Eq.  (3. 23)  it  follows  that  Q*'S  is  antisymmetric. 
Thus  since  Q*1  is  symmetric  the  impedance  Z,  Eq. 
(3.18),  is  Hermitian, 

Z  =  (Zr)*,  (3.26) 

where  *  indicates  complex  conjugation. 

IV.  DEGENERACY 


The  considerations  in  Sec.  Ill  were  largely  based  on 
the  fact  that  the  form  a  complete  set  when  all  pa  are 
different.  However,  degeneracy  may  occur.  For  de¬ 
finiteness,  suppose  that 

p\*p2  =/>3 

,  p<=/>r.  p5=/>?.  (4.D 

Pi*l>S=h 

The  degeneracy  may  either  occur  for  reasons  of  sym¬ 
metry  or  accidentally  for  a  definite  velocity. 

If  the  degenerate  N  still  possesses  six  independent 
eigenvectors,  the  considerations  of  Sec.  Ill  apply  in  toto. 

However,  with  degeneracy  one  is  not  assured  that  six 
independent  eigenvectors  exist  as  in  the  nondegenerate 
case,  for  which  theorem  (3.4)  applies.  When  the  de¬ 
generate  N  does  not  possess  six  independent  eigenvec¬ 
tors,  we  say  that  we  have  the  case  of  nonsimple  de¬ 
generacy.  Static  deformation  [i.e. ,  v  =  0  in  Eq.  (2.6)] 
in  elastically  isotropic  bodies  is  such  a  case.8 

In  the  case  of  nonsimple  degeneracy,  one  forms  a 
complete  set  by  the  introduction  of  generalized  eigen¬ 
vectors  associated  with  the  degenerate  eigenvectors. 
Suppose  the  degenerate  N  corresponding  to  the  example 
(4. 1)  only  allows  four  independent  eigenvectors  £,,  £z, 

4,,  and  £s.  Generalized  eigenvectors  £3  and  as¬ 
sociated  with  the  degenerate  eigenvectors  £2  and  can 
then  be  formed  according  to  the  scheme 


will  be  called  the  surface  impedance.  A  similar  concept 
was  introduced  by  Ingebrigtsen  and  Tonning, 3  but  in  an 
eigenvector  representation.  In  Sec.  IV,  the  above  in¬ 
tegral  expression  for  Z  will  be  shown  to  be  valid  also 
in  cases  of  degeneracy. 

The  matrices  S,  Q,  and  B  introduced  in  the  foregoing, 
are  not  independent,  but  obey  certain  completeness 
relations.  From  Eq.  (3. 12),  applying  N  twice  gives 

N2«a=-?0,  (3.19) 

and  since  the  £a  form  a  complete  set 

(3.20) 

where  {  is  an  arbitrary  six-dimensional  vector.  Thus 
^+(qJ)=°,  (3.21) 


(N-/>.,I)$3=$2,  (4.2) 

(N-/>sI)?g=?5.  (4.3) 

General  theorems  exist5’6  stating  tliat 

the  eigenvectors  together  with  their  associated  ^ 
generalized  eigenvectors  form  a  complete  set. 

Let  |s(y)  be  the  condition  at  depth  y  below  the  surface 
when  the  force-displacement  condition  at  the  surface  is 
£3,  and  similarly  let  £2(y)  correspond  to  £2  at  the  sur¬ 
face.  Combination  of  Eqs,  (4.2)  and  (2.8)  show,  since 
N  and  N-/>2I  commute,  that 

(N-/>2I)«3(y)  =  «2(y).  (4.5) 

£3(y)  remains  the  generalized  eigenvector  associated 
with  |2(y).  The  explicit  y  dependence  of  $3(y)  is  obtained 
by  integration  of 
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gp  ^)=imiiy)=ik{j>AM  +  hfy)\,  (4.6) 

which,  since 

{2(y)  =  exp tikp#)Zt,  (4.7) 

yields 

.  ijiy)  =  exp(i7i/>jy)  $3  +  iky  exp()A’/)2y)  (2.  (4.8) 

It  is  apparent  that  since  £2(y)  is  decaying  with  in¬ 
creasing  y  because  lmp2  >0,  ?3(y)  will  also  vanish  as 
y  —  “.  Thus  (j3(y)  is  also  an  allowed  component  in  a 
surface-wave  construction. 

Next  consider  £3  as  a  function  of  the  angle  in  N(</>), 
£3(</>).  {3  is  not  invariant  to  rotations  0  and  the  question 
arises  whether  Eq.  (3.12),  which  is  the  basis  for  the 
integral  impedance  expression,  holds  for  {3. 

f3(4>)  obeys  the  equation 

[N(«)-/>j(0)I]«3(0)  =  {j,  (4.9) 

where  $2  is  independent  of  <p.  Differentiating  this  equa¬ 
tion  with  respect  to  <j>  and  using  Eqs.  (2. 13)  and  (3.6) 
(which  holds  for  p2),  we  obtain 


(NJ-/>22I)«3  =  (N-/>2I)^-  «„ 


(4. 10) 


from  which  follows,  by  combination  witli  Eq.  (4.9), 

(N-/>2D^«3  =  2/>j«2.  (4.11 

Comparing  Eqs.  (4.11)  and  (4.9)  one  sees  that 

Jq  ?3  =  2/>2£3  +  y£2.  (4.12 

where  y  is  constant.  Actually  £3  is  not  unique;  it  is 
unique  on  to  within  an  arbitrary  component  of  {2.  We 
may  choose  y  =  0.  Integration  of  Eq.  (4. 12)  then  yields 


{3(<#,)  =  ?3(°)exp[K(^)], 


where 


(4. 14) 


k(<P)  =  2  J*  p2(<t>)d<p. 

Consider  now  N(0)$3(O).  By  Eqs.  (4. 9)  and  (4. 13) 
N(*){,(0)  =  exp[-  kW) Nfo)  «,(«) 

=  exp[—  k(<#»)] [/>2(^»)  £3(0)  + 12] 

-PtW  &>(°)  +  exp[-  k(4>)U2. 

Thus  making  angular  averages 
‘  N{3(0)=i?3(0)  +  exp[-  x(</))]£2. 

From  Eqs.  (4. 15)  and  (3.6) 


From  Eq.  (3.6) 

I7pi=o, 

so  finally 

exp[-  k(<I>)\  =  0. 


(4. 16) 


(4. 17) 


(4.20) 


exp[-  «(</)))  =  0.  (4.21) 

[The  special  case  p2(0 )  =  »'  must  be  considered  separate¬ 
ly,  but  the  conclusion  is  the  same.  ] 

Equation  (4. 17)  thus  reduces  to 

NU0)  =  /y0)  (4.22) 

and  we  have  readied  the  important  conclusion  that  Eq. 

(3. 12)  applies  both  to  eigenvectors  and  associated  gen¬ 
eralized  eigenvectors.  Thus,  Eq.  (3. 14)  also  applies  to 
eigenvectors  as  well  as  generalized  eigenvectors.  Re¬ 
calling  theorem  (4.4)  as  the  equivalent  of  theorem  (3.4), 
it  also  follows,  by  reasoning  identical  to  that  of  Sec.  Ill, 
that  A|,  A2,  and  A3  form  a  complete  set  in  the  three- 
dimensional  space  and  that  Z,  Eq.  (3. 18),  still  deter¬ 
mines  what  force  J  at  the  surface  is  needed  to  maintain 
a  given  amplitude  A.  The  integral  expression  (3. 18)  for 
the  surface  impedance  is  generally  valid  for  v<uL. 

It  is  now  also  immediately  clear  that  Eqs.  (3. 19)— 

(3.  26)  remain  valid  with  nonsiniple  degeneracy. 

V.  SOME  PROPERTIES  OF  THE  IMPEDANCE  AND 
SOME  PROPERTIES  OF  THE  LIMITING  WAVE 

We  consider  the  velocity  range  0<v<vL  for  which  the 
theory  in  the  foregoing  applies.  The  appearance  of  the 
limiting  wave  at  v  =  vL  will  also  be  discussed  in  more 
detail. 

The  impedance  Z,  Eq.  (3. 18),  is  Hermitian.  Thus  it 
may  be  represented  in  the  dyadic  form 

Z  =  f}  A,e, .  ef,  efeJ  =  6lj,  (5,1) 

where  e{  are  complex  eigenvectors  obeying  orthogonality 
relations  and  A(  are  real  eigenvalues. 

The  Lagrange  function  per  unit  length  per  unit  depth 
lor  a  surface  disturbance  in  uniform  motion  is 

L  =  'A  J’A  +  JA’),  (5.2) 

or  by  Eqs.  (3.18)  and  (3. 26), 

-L=- [A*ZA.  (5.3) 

At  ii  -  0,  -  L  is  strain  energy  only  and  thus  positive  de¬ 
finite.  Since  A  in  (5.3)  is  arbitrary,  it  follows  that 

Aj  >0,  is  1,2, 3  when  v-Q.  (5.4) 

Keeping  the  amplitude  A  fixed  but  changing  the  ve¬ 
locity,  the  relation 


—m, 


is  valid,  where  T  is  kinetic  energy  per  unit  length  per 
unit  depth.  Thus,  by  Eq.  (5. 3) 


*(</>)  =  - 2  f 


pdp  _  .  1  +/>?(</>) 
l+p  l+/4(0) 


W •  3l,A'(l;z)A=‘2r' 


(4.19) 


Since  T  is  positive  definite,  it  follows  that 


—  Z  is  negative  definite 


and 

X{,  t  =  1 , 2, 3  decrease  with  increasing  i>,  (5.8) 

A  similar  conclusion  based  on  the  relation 


where  E  is  the  total  energy,  was  reached  by  Ingebngtscn 


and  Tonning. 3 

From  Eqs.  (3.18)  and  (3,23)  it  follows  that 

AZA  =  -ArAQ*'A.  (5.9) 

Since  Q‘‘  is  real  and  negative  definite  when  t  "vL,  it 
follows  that 

AZA>0  for  real  A,  (5.10) 

the  possible  equality  in  (5. 10)  referring  to  which 

will  be  discussed  later. 

The  power  flux  per  unit  area  in  they  direction,  VFy,  is 

Wy=(u>/4i)(J*A- JA*).  (5.11) 

Of  course,  in  a  surface  disturbance  in  uniform  motion 
parallel  with  the  surface 

WV-0.  (5-12) 

The  impedance  (3. 18)  satisfies  this  requirement 

W'y=  -(<u/4i){A*[(Zr)*-ZlA}  =  0.  (5.13) 


imaginary,  it  now  follows  from  Eq.  (5.11) 

Jt  At  =  0.  (5.18) 

VI.  EXISTENCE  THEOREMS 

Most  of  the  facts  listed  in  Sec.  ire  well  known  and 
therefore  detailed  proofs  have  been  omitted.  The  only 
property  which  appears  not  to  be  commonly  realized  is 
that  of  (5. 10),  which  also  will  be  of  some  importance 
in  the  considerations  to.follow. 

However,  what  apparently  has  not  been  realized  be¬ 
fore  is  that  the  preceding  list  of  facts  provide  sufficient 
basis  for  the  construction  of  existence  theorems. 

If  a  velocity  vR<v L  for  which 

IIZ(«*)||=0  '  (6.1) 

can  be  found,  then  a  surface  wave  satisfying  the  con¬ 
dition  of  free  surface,  i.e. ,  a  Rayleigh  wave  (RW),  with 
phase  velocity  vR,  exists.  Condition  (6. 1)  is  the  con¬ 
dition  that  J  =  0  in  Eq.  (3. 17)  be  compatible  with  some 
nonzero  A.  The  amplitude  of  the  surface  wave  satisfies 

Z(n«)A/(=:0.  (6.2) 

Condition  (6. 1)  is  equivalent  with  the  condition  that 
one  of  the  eigenvalues  X,  in  Eq.  (5. 1)  be  zero  at  v  =  vR. 

Comparing  (5. 10)  and  Eq.  (6.2),  it  is  immediately 
clear  that  the  RW,  if  it  exists,  is  elliptically  polarized. 


As  v  —  Vt,  Q  diverges  in  the  sense  that  one  of  its 
eigenvalues  approaches  infinity.  Q"1  is  convergent,  and 
one  of  the  eigenvalues  approaches  zero.  S  diverges, 
but  Q'‘S  converges.  This  is  ensured  by  the  eigenvalue 
in  Q'1  approaching  zero.  Thus  as  is  well  established 
and  can  be  ascertained  in  the  present  formalism  by  the 
methods  of  Barnett  and  Lothe, 1  Z (/»)  approaches  a  finite 
limit  as  v  —  vL,  and  we  define  7,(vL)  through  this  limit 

Z(vi)  =  lim  Z(»>).  (5.14) 

V~VL 

The  limiting  wave  arises  by  coalescence  of  two  partial 
waves  ?a  and  as  v  —  vL.  Since  Z{v)  applies  to  one  of 
the  two  coalescing  waves  (that  for  which  lm/>a  >0),  we 
must  expect  Z(vL)  to  apply  to  the  limiting  wave  itself 

Jz  =  -Z(yJAx,.  (5.15) 


Now  make  the  following  supposition:  The  limiting  wave 
does  not  satisfy  the  condition  of  free  boundary 

Ji^O  (6.3) 

and  no  RW  exists  for  v<vL.  We  shall  see  that  a.con- 
tradiction  arises. 

From  Eq.  (6.3) 

J£Jl>0,  (6.4) 


or,  by  Eq.  (5.15),  (5.16),  and  Eq.  (3.26) 
Aj,Zz((;t)Ai>0, 

which  by  Eq.  (5. 1)  takes  the  form 

1*1 


(6.5) 

(6.6) 


Indeed,  the  validity  of  this  relation  can  be  proven. 

A i  is  the  eigenvector  of  a  real  matrix,  Eq.  (2. 15). 
Thus  kL,  the  amplitude  of  the  limiting  wave,  can  be 
chosen  as  purely  real .  The  limiting  wave  is  linearly 
polarized.  We  will  always  choose 

Ax,  real.  (5.16) 

With  At  real  and  the  fact  that  the  limiting  wave  is  a 
bulk  wave,  it  follows  that 

Jx,  purely  imaginary.  (5. 17) 

From  the  construction  (see  for  example  Ref.  1)  of 
the  limiting  wave  from  the  slowness  surface,  it  follows 
that  the  group  velocity  is  parallel  with  the  surface  so 
that  Wy  =  0  also  for  the  limiting  wave.  Alternatively,  we 
may  say  that  the  limiting  wave  has  inherited  this  proper¬ 
ty  from  Im/>a>0.  Since  kL  is  real  and  3L  purely 


where  A,^  is  the  value  of  X,  at  v  =  vc,  etc. 

According,  to  Eq.  (5. 18) 

£/AliJell*Aj,  =  0.  (6.7) 

»■! 

Let  AP)x,  be  the  largest  eigenvalue.  Multiply  Eq.  (6.7) 
by  X Pti  and  subtract  it  from  inequality  (6.6) 


Ax,£,(X4  £,  —  \p  L)  |e,Ax, 
f»i 


!2 


0. 


(6.8) 


With  no  RW  for  v<  t>L,  all  X(t£. 5^  0  according  to  Eqs. 
(5.4),  (5.8),  and  (6.1).  All  terms  X(ii-  Xyii  are  negative 
or  zero  according  to  the  definition  oi  Inequality 
(6.8)  is  therefore  incompatible  with  our  assumption  and 
we  have  arrived  at  the  following  theorem: 

When  the  limiting  wave  does  not  satisfy  the  con- 
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(7.4) 


dition  of  a  free  surface,  it  can  be  asserted  that 
a  Rayleigh- wave  solution  exists  for  some  vR<vL. 

(6.9) 

Suppose  two  RW  solutions  exist,  for  vRi<vL  and 
vxi<vJ.7  respectively.  Then  two  eigenvalues,  say  X,  and 
X2,  must  be  negative  at  vL, 


X2fX,<0. 

We  can  always  find  a  real  vector  A  which  is  perpendicu¬ 
lar  to  both  the  real  and  imaginary  part  of  e3li.  For  this 
real  vector  A 


^2=  (C; ;o/p)l/2, 

V,=  (Cn/p)u\ 

where  p  is  the  density.  Only  the  wave  with  velocity  Fj 
satisfies  the  condition  of  free  surface.  We  are  interested 
in  the  case  when  this  wave  defines  the  limiting  velocity, 
thus 

«-5> 

Since  the  only  restriction  on  C44  is  that  C44  must  be 
positive,  this  is  a  possible  case. 


AZA— XJtx. I ei,x.A |2  X2, x, | e2, /.A  |2 < 0,  (6.11) 

in  violation  of  inequality  (5. 10). 

Thus  we  also  have  the  theorem: 

There  is  at  most  one  RW  solution  for  v<vt.  (6. 12) 

VII.  DISCUSSION 

Some  complications  arise  when  more  than  one  limiting 
wave  exist  at  vL.  The  various  cases  that  may  arise  have 
previously  been  discussed  quite  exhaustively  by  the 
present  authors.  The  general  conclusion  is  a  theorem 
actually  somewhat  less  restrictive  than  theorem  (0.9): 

When  one  of  the  limiting  waves  at  vL  does  not 
satisfy  the  condition  of  a  free  surface,  it  can 
be  asserted  that  a  Rayleigh-wave  solution  exists 
for  some  vH<vL.  (7.1) 

When  the  limiting  wave  does  satsify  the  condition  of 
a  free  surface,  a  Rayleigh  wave  may  or  may  not  occur 
for  some  vK<vL,  In  this  situation  both  cases  occur.  The 
previous  considerations 1  by  the  present  authors  for 
transverse  isotropy  contained  some  algebraic  errors 
and  we  will  use  this  opportunity  to  present  a  corrected 
discussion. 

•We  consider  wave  propagation  in  the  x  direction  on  the 
surface  perpendicular  to  tire  y  axis.  The  elastic  coef¬ 
ficients  Ctl,  C33,  C44,  CG0,  C, 2,  and  Crj  describe  the 
crystal.  Isotropy  about  the  2  axis"  requires  that 

2C66  =  Cu-  C,2.  (7.2) 

Lattice  stability  requires  that  Cm  C44,  and  C00  be 
positive  and  that 

C\\-C\i>Q.  (7.3) 

Three  bulk  waves  enter  the  problem.  Their  respective 
velocities  are 

^1  —  (Pu/p)l,i> 


From  the  isotropy  condition  (7.2)  and  the  lattice 
stability  requiremc.it  (7.3)  it  follows  that 

F2<F3.  (7.6) 

Thus  we  are  discussing  the  case 

F,  <  F2  <  F3.  (7.7) 

Explicit  calculations  show  that  vR  is  a  function  of  F2  and 
F3  only  and  because  of  inequality  (7. 6)  vR  is  smaller 
than  F2,  so  that 

t>*=t'n(F2,F3)<  F2.  (7.8) 

Now  given  F2  and  F3  and  thus  also  vR,  we  may  vary  F, 
independently  by  changing  C44  without  violating  lattice 
stability,  to  realize  either  one  of  the  cases 

>,R<vi  —  V’i  <  F2  <  F3,  (7.9) 

or 

t>L  =  F,  <  i>*  <  F2  <  F3.  (7.10) 

In  case  (7. 10)  there  is  no  RW  for  v  <  vL.  In  this  case 
the  surface  wave  with  velocity  vR  belongs  to  a  branch  of 
leaky  surface  waves. 9 

The  theoretical  conclusions  of  this  paper  are  in  com¬ 
plete  accord  with  and  in  fact  rationalize  the  numerical 
work  reported  by  Farnell. 9 
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